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PREFACE 

The  material  described  in  this  report  is  the  result  of  collaborations, 
over  different  periods  of  time,  of  several  scientists.  Except  as  noted,  all 
are  members  of  the  Health  Physics  Division  of  the  Oak  Ridqe  National 
Laboratory. 

Over  the  course  of  this  project,  the  following  have  taken  leading  roles 
in  guiding  and  performing  this  work: 

E.  T.  Arakawa 
R.  D.  Birkhoff 
R.  H.  Ritchie 

Others  making  substantial  contributions  to  this  project: 

V.  E.  Anderson  (Computer  Sciences  Division,  ORNL) 

J.  C.  Ashley 

A.  J.  Braundmeier  (Consultant  from  Southern  Illinois 
Universi ty) 

T.  A.  Call co tt  (Consultant  from  The  University  of  Tennessee) 

L.  C.  Emerson  (now  Metals  and  Ceramics  Division,  ORNL) 

R.  N.  Hamm 

T.  Inagaki  (now  Osaka  Kyoiku  University,  Japan) 

C.  J.  Tung  (Postdoctoral  Fellow  from  The  University  of  Tennessee) 

M.  W.  Williams 

We  point  out  that  the  Ph.D.  research  of  Dr.  C.  J.  Tung  (Ph.D.,  The 
University  of  Tennessee,  1975)  was  largely  concerned  with  problems  arising  in 
the  course  of  this  project. 

The  continuing  advice  and  assistance  of  E.  A.  Burke,  J.  C.  Garth,  and 
J.  N.  Bradford  of  the  Rome  Air  Development  Center  is  gratefully  acknowledged. 
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I.  PROGRAM  OBJECTIVE 


The  objective  of  this  program  was  to  gain  better  understanding  of  the 
basic  interactions  of  ionizing  radiations  with  materials  of  importance  as 
electronic  components  of  communications  or  guidance  systems.  The  knowledge 
gained  should  lead  to  a better  design  of  components  and  an  appreciation  of 
their  response  to  abnormal  situations  such  as  could  exist  during  the  passage 
of  an  extremely  high  energy  cosmic  ray  or  might  occur  in  the  high  radiation 
environment  associated  with  a nuclear  explosion. 

Our  approach  has  been  both  theoretical  and  experimental.  The  energy 
spectrum  of  electrons  generated  by  ionizing  radiations  and  the  attenuation 
length  or  range  of  electrons  as  a function  of  electron  energy  have  been 
studied.  A spherical,  electrostatic  electron  spectrometer,  the  "Keplertron" , 
specifically  designed  to  measure  slowing-down  spectra  in  solids,  was  used  to 
obtain  the  electron  energy  snectra.  These  spectra  were  then  analyzed  in  terms 
of  existing  theories  and  in  terms  of  new  models  developed  for  this  project. 
Attempts  have  been  made  to  obtain  values  for  electron  attenuation  lengths 
from  photoelectric  yield  observations.  This  has  proved  to  be  particularly 
difficult  for  small  electron  attenuation  lengths  of  10  A:  a conclusion  sup- 

ported by  theoretical  studies  of  photoemission. 

In  line  with  our  overall  objective,  theoretical  models  and  subsequent 
calculations  based  on  experimental  results  have  resulted  in  tabulations  of 
mean  free  path,  energy  loss,  energy  straggling,  range,  and  range  straggling 
for  electrons  in  various  metals,  semiconductors,  and  insulators. 
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II.  REVIEW  OF  PROGRAM 


A.  ELECTRON  SLOWING-DOWN  SPECTRA 

When  a medium  is  irradiated  by  ionizing  radiations,  a flux  of  electrons 
is  generated.  The  primary  flux  may  consist  of  high  energy  Compton  or  photo- 
electrons from  X-ray  irradiation  or  beta  rays  from  radioactive  decay 
processes.  As  these  high  energy  electrons  slow  down  in  the  medium,  they 
excite  other  electrons.  The  more  energetic  of  these  are  called  delta  rays  , 
while  those  with  a few  Rydbergs  or  less  of  energy  are  usually  referred  to  as 
secondaries  . If  the  medium  is  very  large,  homogeneous,  and  isotropic, 
attention  may  be  focused  on  the  generation  and  slowing-down  of  these  electrons 
without  regard  for  the  spatial  variation  of  the  flux. 

The  purpose  of  this  research  was  to  determine  the  spectrum  of  electrons 

which  are  generated  by  ionizing  radiations  in  materials  of  importance  in 

electronics  applications.  Such  materials  include  metals,  semiconductors,  and 

1 8 

insulators.  The  electron  spectra  of  interest  have  been  designated  as 
slowing-down  spectra  or  degradation  spectra  by  theorists,  Spencer  and 
Fano  in  the  first  case,  and  Platzman  in  the  second.  Our  study  has  uti- 
lized an  experimental  facility  at  the  Oak  Ridge  National  Laboratory  called 

1 9 

the  "Keplertron" , which  is  a spherical,  electrostatic,  electron  spectrometer 
and  sources  made  radioactive  in  the  Oak  Ridge  Research  Reactor  or  the  High 
Flux  Isotope  Reactor.  A companion  theoretical  effort  has  involved  numerical 
solutions  of  the  Boltzmann  transport  equation  and  comparisons  with  the 
experimental  results. 

The  materials  studied  in  this  portion  of  the  contract  were  the  semicon- 
ductor Si  and  the  insulator  Al^O^.  Investigations  of  these  materials  formed 
a natural  extension  of  earlier,  comprehensive  studies  carried  out  in  our 
group  on  several  metallic  media.  Electron  slowing-down  studies  on  metals 
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have  been  reviewed  in  a publication  of  the  National  Academy  of  Sciences  and 
elsewhere. ' 

For  the  semiconductor  Si,  the  electron  slowing-down  flux  spectrum  has 

been  measured  over  the  energy  range  from  a few  eV  to  30  keV  using  neutron- 
2 

activated  samples.  The  spectrum  is  found  to  be  similar  to  slowing-down 
spectra  generated  in  metals  by  beta  particles.  Absolute  agreement  with 
predictions  of  the  Spencer-Fano  theory1^  is  found  for  energies  above  10  keV. 
Experimental  fluxes  are  found  to  be  somewhat  larger  than  these  predictions 
for  lower  energies.  A companion  theoretical  calculation  has  been  carried  out 
based  on  valence  electron  cross  sections  obtained  from  an  electron  gas  model 
and  on  inner-shell  cross  sections  obtained  from  sum-rule-corrected  binary 
encounter  theory.*  These  cross  sections  describe  the  effects  of  plasmon 
creation  and  electron-hole  pair  production  in  the  valence  band  and  excitation 
of  electrons  from  the  K-  and  L-levels.  Monte  Carlo  calculations  of  the 
slowing-down  spectrum  were  made  using  these  cross  sections.  This  theory,  as 
well  as  an  empirically  she! 1 -corrected  Spencer-Attix  theory,  gives  a spec- 
trum with  magnitudes  which  are  smaller  than  experimental  values  by  a factor 

of  % 4 in  the  10  to  100  eV  range.  A detailed  account  of  this  work  has  been 

3 

published  in  The  Physical  Review. 

Since  insulators  are  frequently  used  in  components  of  electronic  devices, 
it  is  important  to  know  their  response  to  radiation.  Aluminum  oxide  was 
chosen  as  a representative  material  for  this  study  partly  because  of  its  use 
in  M0S  devices  where  it  replaces  silicon  dioxide  in  those  casv:'  where 


*A  description  of  our  theoretical  work  on  classical  binary  collision 
cross  sections  of  atomic  systems,  which  has  been  prepared  for  publication,  is 
included  in  this  report  as  Appendix  A. 
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radiation  hardening  is  important.  Details  of  sample  preparation  and  proce- 
dures for  the  experimental  measurements  are  described  in  Appendices  B-l  and 

164 

B-2  along  with  the  observed  electron  slowing-down  spectrum  for  Dy  beta 
rays  in  Al^O-j.  The  measured  spectrum  exhibits  an  overall  shape  which  is 
similar  to  that  for  the  semiconductor  Si  and  the  metals  measured  previously. 
The  electron  flux  is  generally  higher  than  that  seen  for  a typical  metal. 

Two  theoretical  calculations  have  been  carried  out  to  compare  with  the 

3 

measured  fluxes.  One  calculation,  analogous  to  that  performed  for  silicon, 
used  Monte  Carlo  techniques  employing  inelastic  cross  sections  obtained  from 
a sum-rule-constrained,  classical  binary  collision  model.  The  resulting 
spectrum,  shown  in  Appendix  B-2,  has  been  described  previously.  ’ Above 
10  keV  the  calculated  spectrum  shows  approximately  the  same  energy  dependence 
as  the  data  but  is  % 40  lower.  At  lower  energies  the  theory  falls  somewhat 
below  the  experimental  data  but  tends  to  converge  with  it  at  m 10  eV.  The 
overall  agreement  was  encouraging  considering  the  uncertainties  involved  in 
using  the  approximate  inner-shell  cross  sections  over  such  a broad  range  of 
energies. 

In  hopes  of  resolving  the  rather  large  discrepancy  between  theory  and 
experiment  in  the  energy  range  from  around  40  eV  to  2000  eV,  a second  theore- 
tical calculation  was  done.  A theoretical  model  for  the  response  function  of 
an  insulator  was  developed  (see  Appendix  C)  to  describe  the  behavior  of  the 
valence  band  electrons  and  then  applied  to  A 1 , 0 ^ . Cross  sections  for  excita- 
tion of  inner-shell  electrons  from  the  A1  and  0 ion  cores  based  on  detailed, 
theoretical  calculations  of  atomic  generalized  oscillator  strengths,  combined 
with  a tight  binding-OPW  model  description  of  valence  band  electrons,  were 
used  in  a Bol tzmann-type  integral  equation  to  obtain  the  electron  slowing- 
down  flux.  These  results  are  described  briefly  in  Reference  6.  Later  this 
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calculation  was  further  improved  by  including  Auger  cascade  transitions. 

This  more  detailed  calculation,  summarized  in  Appendix  D,  gave  essentially 
the  same  result  for  the  electron  slowing-down  spectrum  in  Al^O-  as  the 
earlier  Monte  Carlo  calculation  and  thus  did  not  resolve  the  discrepancy 
noted  above  between  theory  and  experiment. 

A recent  review  paper  on  electron  slowing-down  spectra5  gives  an  over- 
view of  the  studies  on  Si  and  Al^O^  described  above  as  well  as  the  studies  of 
metallic  media.  A detailed  account  of  the  theoretical  work  involved  in  the 
electron  slowing-down  flux  calculations  has  been  published  as  an  ORNL  report.7 

B.  ELECTRON  ATTENUATION  LENGTHS  FROM  PHOTOELECTRIC  YIELD  MEASUREMENTS 

A knowledge  of  the  range  of  electrons  as  a function  of  their  energy  in 
each  of  the  materials  involved  in  a communications  system  is  necessary  in 
order  to  estimate  the  effects  of  ionizing  radiations  on  that  system.  Since 
the  range  and  mean  free  path  of  electrons  in  a medium  are  related  to  the 
attenuation  length  associated  with  electrons  emitted  from  the  surface  of  the 
medium,  measurements  of  electron  attenuation  lengths  can  yield  the  desired 
values  of  the  range.  Conventionally  these  quantities  have  been  determined 
from  electron  transmission  measurements  through  thin  films.  However,  for 
electrons  with  energies  of  only  a few  electron  volts  above  the  Fermi  level  up 
to  n,  1000  eV,  it  is  known  that  electron  attenuation  is  very  strong  and  the 

o 

associated  lengths  are  therefore  very  short;  in  fact  £ 10  A.  It  is  thus 
impossible  at  present  to  prepare  homogeneous  films  which  are  sufficiently 
thin  for  these  measurements. 


We  had  previously  proposed^3  a method  of  obtaining  short  electron  attenu- 
ation lengths  from  measurements  of  photoelectric  yield  with  an  analysis  based 

22 

on  a theory  by  Pepper.  Under  this  contract  we  used  our  proposed  method 

to  obtain  values  of  electron  attenuation  lengths  of  photoexcited  electrons 
in  thick  films  of  carbon8  (see  also  Appendix  E-l), 
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of  aluminum,  and  of  polystyrene  from  measurements  of  total  photo- 

electric yield  as  a function  of  photon  angle  of  incidence,  over  the  range  of 
photon  energies  from  20  to  64  eV.  For  a given  photon  energy  the  yield,  Y (e), 
for  a photon  angle  of  incidence  o , relative  to  the  yield  Y ( 0 ) for  normal 
incidence,  shows  a single  maximum  when  plotted  vs  e if  the  film  is  infinitely 
thick  and  if  its  refractive  index,  n,  is  less  than  unity.  The  position  of  this 
maximum  in  G(e),  where  G ( e ) = Y(e)/Y(0),  is  in  the  vicinity  of  ec  where  ec  = 
sin  ^n,  when  n is  close  to  unity.  The  magnitude  of  G(ec)  should  be  sensitive 
to  both  the  extinction  coefficient,  k,  and  the  electron  attenuation  length, 

1.  Initially,  a three  parameter  least-squares  fit  of  the  photoelectric  yield 
to  the  theory4  was  attempted.  However,  it  was  soon  apparent  that 

we  could  not  obtain  meaningful  values  of  both  k and  L from  a least-squares 
fit  of  an  experimental  curve  of  G(o)  vs  6.  Values  of  k were  thus  obtained 
independently  and  a two  parameter  least-squares  fit  of  G(o)  vs  e employed  to 
yield  values  of  n and  L.  For  carbon,  thin  self-supporting  films  could  not  be 
prepared  for  transmission  measurements  and  k values  were  obtained  from  obser- 
vations of  reflectance  as  a function  of  angle  of  incidence  using  thick 
evaporated  carbon  films  similar  to  those  used  for  the  photoelectric  yield 
measurements.  The  uncertainties  in  the  k values  obtained  by  this  technique 
were  quite  large  over  the  range  of  photon  energies  from  20  to  64  eV.  However, 
the  reflectance  vs  photon  angle  of  incidence  measurements  yield  both  n and  k 
values  and  it  was  felt,  at  that  time,  that  a consistency  check  could  be  made 
by  comparing  the  values  of  n obtained  at  a given  photon  energy  from  the 
reflectance  and  from  the  photoelectric  yield  measurements.  We  later  found 
that  the  consistency  we  obtained  in  n for  carbon  from  the  two  experimental 
techniques  was  not  a sufficient  condition  for  reliability  of  the  L values 
obtained.  For  polystyrene,  thin  films  suitable  for  optical  transmission 
measurements  over  the  whole  energy  range  could  be  fabricated.  The  values  of 
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k obtained  have  formed  the  basis  for  a detailed  study  and  analysis  of  the 
optical  properties  of  polystyrene^  (see  Appendix  F).  These  optical  proper- 
ties are  required  in  the  calculations  involved  in  producing  a tabulation  of 
electron  energy  loss,  energy  straggling,  range,  and  range  straggling  for 
electrons  in  polystyrene  (compare  Part  II-C  and  Part  III).  Despite  the  con- 
sistency of  the  n values  obtained  from  reflectance  and  from  photoelectric 
yield  data,  the  magnitude  of  some  of  the  L values  obtained  for  carbon  were 
considered  questionable,  and  a detailed  sensitivity  analysis  was  thus  under- 
taken. Three  techniques  had  been  used  to  obtain  k values  for  carbon: 
transmission  through  thin  carbon  films  on  a thin  A1  substrate,  reflectance  as 
a function  of  photon  angle  of  incidence  using  thick  carbon  films,  and  the 
three  parameter  least-squares  fit  to  the  G ( e ) vs  9 curves.  The  values  used 

O 

in  the  calculations  of  L were  those  obtained  from  the  reflectance 
measurements.  These  were  thought  to  be  the  most  accurate  and  were  in  fact 
the  average  of  all  the  results.  Uncertainties  in  k of  % 50%  were  estimated 
from  the  range  of  k values  obtained  by  the  three  methods.  On  the  other  hand 
G(e)  was  found  to  be  reproducible  to  'v  1%.  These  uncertainties  in  k and 
G(e)  were  found  to  lead  to  large  uncertainties  AL(k)  and  aL(G),  respectively, 
in  the  values  calculated  for  L from  the  experimental  G(o)  vs  9 data.  For 
example  at  22.3  eV,  L = (113  + 91)  A where  AL(k)  = 81  A and  AL(G)  = 10  A,  at 
40.8  eV,  L = (19  ± 30)  A where  AL( k)  = 29  A and  aL(G)  = 1 A,  and  at  59.2  eV, 

L = (10  ± 45)  A where  AL(k)  = 44  A and  AL(G)  = 1 A.  Appropriate  error  bars 
should  be  added  to  the  L values  graphed  for  carbon  in  Reference  8.  If 

o 

L is  small  (£  10  A),  as  is  generally  thought  over  the  whole  of  the  energy 
range  studied  here  for  carbon,  it  can  be  shown  that,  in  general,  to  obtain 

O 

L to  +1  A (i.e.,  to  at  best  ^ +10%),  not  only  G(9)  but  also  k must  be  known 
to  better  than  +V  . There  may  be  other  inaccuracies  in  the  values  of  L 
calculated  by  this  technique  due  to  the  sample  not  being  homogeneous  and/or 
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the  sample  surface  not  being  perfectly  clean,  smooth,  and  planar.  In  addition 
the  model  used  in  the  derivation  of  the  theoretical  variation  of  6(0)  with  0 
in  terms  of  the  definition  of  L may  be  inadequate.  In  fact  photoelectric 
yield  continues  to  be  the  subject  of  theoretical  and  experimental  investiga- 
tions designed  to  elucidate  the  processes  involved  and  the  significance  of 
the  measured  electron  attenuation  length. 

Since  analysis  of  6(0)  vs  0 yielded  questionably  high  values  of  L (%  150  A 

o 

to  50  A)  for  carbon  in  the  photon  energy  region  from  20  to  30  eV,  we  turned 
to  a second  technique11  (see  Appendices  E-2  and  G)  based  on  photoelectric 
yield  measurements.  A wedge-shaped  carbon  film  was  deposited  on  a glass  sub- 
strate and  the  thickness,  t,  as  a function  of  position  along  the  wedge 
obtained  from  a combination  of  transmi ttance  and  independent  interferometric 
measurements.  The  total  yield  Y(t)  normalized  to  the  yield  Y(°°)  for  a thick 
sample,  was  obtained  as  a function  of  t over  the  same  energy  range  of  20  to 
64  eV.  These  data  indicated  that  the  values  of  L for  photoexcited  electrons 
in  carbon  are  much  smaller  in  the  region  of  20  to  40  eV  than  those  obtained 
from  observations  of  G(e)  vs  0.  All  electron  attenuation  lengths  in  solid 

o 

media  are  thought  to  be  small  (£  10  A)  in  the  region  of  % 15  to  ^ 1000  eV. 

The  values  inferred  by  inspection  of  Y(t)/Y(<»)  vs  t were  in  fact  in  this 

range  for  carbon.  However,  a detailed  interpretation  of  these  data  for 

evaporated  carbon  films  was  subject  to  several  theoretical  and  experimental 

uncertainties.  Gold  is  simpler  to  work  with  and  the  optical  properties  of 

gold  are  better  documented  in  the  literature  than  those  of  carbon.  Hence,  it 

was  thought  that  measurements  of  Y(t)/Y(™)  vs  t for  gold  films  would  provide 

a check  on  the  carbon  data  and  its  interpretation.  Data  was  thus  obtained 

22 

on  gold  films  and  an  analysis  in  terms  of  the  theory  proposed  by  Pepper 
was  attempted.  The  range  of  film  thickness  required  for  an  accurate  determi- 
nation of  attenuation  length  depends  on  the  optical  properties  of  the  film 
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| and  the  value  of  L of  the  emitted  electrons.  For  a strongly  absorbing  film 

and  small  L the  thickness  available  should  range  from  t < L to  t = <*>.  For 
weakly  absorbing  films,  one  does  not  need  to  have  films  with  t < L,  but 
small  enough  thicknesses  should  be  used  that  the  first,  and  strongest,  struc- 
ture in  Y(t)/Y(°°)  vs  t is  accurately  specified.  Furthermore,  these  films 
need  to  be  homogeneous  and  continuous  and  to  have  smooth  interfaces  with  the 
substrate  and  with  vacuum  in  order  for  Y(t)/Y(°°)  vs  t to  yield  meaningful  L 
values.  A literature  survey  has  now  shown  that  using  currently  available 
techniques  for  making  thin  films,  all  films  are  discontinuous  for  t less 

o 

than  about  50  to  300  A,  the  thickness  at  which  continuity  occurs  depending 
on  the  method  of  preparation.  In  many  cases,  values  of  L calculated  from 
observations  on  ultrathin  films  would  be  more  dependent  on  the  method  of  film 
preparation,  and  hence  film  structure,  than  on  the  basic  electronic  properties 
of  the  material  of  the  film.  Thus  this  method  is  generally  unsuitable  for 
the  determination  of  L values  in  the  UV  and  soft  X-ray  region.  In  view  of 
this  conclusion,  some  L values  obtained  by  this  method  and  published  in  the 
literature  turn  out  to  be  meaningless,  except  insofar  as  they  set  upper 
bounds  on  the  actual  values  of  L.  This  was  the  case  for  our  Y(t)/Y(=>)  vs  t 
data  for  both  carbon  and  gold. 

Examination  of  variations  of  these  two  methods,  such  as  6(e)  vs  e for 
thin  films  of  known  thickness,  or  the  ratio  of  photoelectric  yields  for  a 
thin  film  illuminated  from  the  front  and  from  the  back  through  the  substrate, 
as  a function  of  film  thickness,  reveals  that  all  methods  for  determining  L 
which  involve  photoelectric  yield  measurements  appear  subject  to  the  same 
limitations  when  L is  small  and  fe  material  is  moderately  to  strongly 
absorbing. 

To  summarize  our  attempts  to  get  electron  attenuation  lengths  from 
photoelectric  yield  measurements;  when  the  medium  is  strongly  absorbing. 
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techniques  involving  thick  films  must  be  used  to  determine  L values.  However, 

O O 

to  obtain  L to  ±1  A when  L is  less  than  about  10  A requires  that  k be  known 
to  within  ±1%.  For  strongly  absorbing  materials  in  the  energy  range  of 
interest,  k cannot  presently  be  measured  to  this  accuracy.  Thus  L cannot  be 
obtained  with  any  meaningful  degree  of  accuracy.  When  the  medium  is  weakly 
absorbing  and  L is  short,  both  thick-  and  thin-film  techniques  can  be  used, 

but  data  should  be  interpreted  with  great  care.  Films  should  not  be  employed 

that  are  so  thin  that  they  are  not  continuous.  Effective  L values  calculated 
from  observations  on  such  films  are  dependent  on  film  structure  and  hence  are 
not  characteristic  of  the  bulk  material.  When  effective  L values  can  be 

obtained  by  a method  employing  thin  films  that  are  thick  enough  to  be  continu- 

ous and  homogeneous,  they  should  be  corrected  for  size  effects  so  that  the 
electron  attenuation  lengths  obtained  correspond  to  those  for  electrons 
escaping  across  a planar  interface  between  the  semi-infinite  medium  and 
vacuum.  Using  currently  developed  experimental  techniques  L values  can  be 
obtained  for  weakly  absorbing  materials,  but  it  is  extremely  difficult  to  get 
values  accurate  to  even  +10%  by  any  knov/n  photoemission  technique.  Finally, 
especially  for  short  attenuation  lengths,  the  theories  of  photoemission 
employed  in  their  derivation  need  to  be  refined. 

Although,  as  already  described,  observations  of  G(o)  vs  9 from  a thick 
sample  cannot  be  analyzed  for  electron  attenuation  lengths  when  the  sample  is 
strongly  absorbing,  the  relative  yields  for  a given  photon  energy  are  found 
to  be  characteristic  of  the  material.  Generally,  G(p)  vs  0 is  reproducible 
to  + VI  and  is  relatively  insensitive  to  surface  conditions.  These  data, 
currently  unavailable  in  the  open  literature,  are  of  practical  interest  in 
the  development  of  improved  photoemitters  in  the  vacuum  ultraviolet  and  soft 
X-ray  regions.  Our  accumulated  G(o)  vs  9 data  on  a variety  of  metals,  semi- 
conductors, and  insulators^  are  currently  being  prepared  for  publication 
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without  attempting  to  obtain  the  associated  L values.  In  some  cases,  where 
the  material  is  strongly  absorbing,  analysis  of  G( e ) vs  9 assuming  a range  of 
small  L values  may  be  the  easiest  and  most  accurate  way  of  obtaining  the 
optical  properties  n and  k.  This  possibility  is  presently  being  investigated. 

C.  TABULATIONS  OF  ELECTRON  INTERACTION  PARAMETERS  IN  SOLIDS 

A quantitative  description  of  the  interaction  of  electrons  with  matter 

over  a large  range  of  energies  is  a subject  of  basic  importance  in  a wide 

variety  of  theoretical  and  applied  areas.  From  the  theoretical  standpoint, 

calculations,  by  other  groups,  of  energy  loss  and  range  of  electrons  in  many 

different  materials  had  previously  formed  the  basis  of  at  least  two  extensive 
23  24 

tabulations.  ’ Both  of  these  works  were  restricted  to  electron  energies 
keV  and  were  based  on  the  Bethe  theory  of  stopping  power  including  various 
modifications  and  corrections  (e.g.,  density-effect  corrections).  We  felt 
that  calculations  of  mean  free  paths,  energy  loss,  etc.,  for  electrons  of 
energy  <J0  keV,  based  on  the  best  currently  available  theoretical  models, 
would  provide  useful  guides  for  interpretation  of  experimental  data  as  well 
as  input  for  calculations  in  applied  areas.  Theoretical  models  developed  and 
tested  in  our  group  and  elsewhere  over  many  years  have  formed  the  basis  for 
our  calculational  efforts  in  this  part  of  the  program.  Our  attention  has 
focused  on  inelastic  interactions  in  which  electron  excitations  play  a major 
role.  Thus  we  consider  incident  electron  energies  from  a few  eV  to  10  keV. 
Inelastic  processes  in  the  "subexcitation"  region  of  energies  (a  few  eV  or 
less),  for  example  phonon  generation,  must  be  treated  using  different 
theoretical  techniques.  For  the  electron  energy  range  studied  in  this  work, 
the  principal  mechanisms  determining  the  interaction  parameters  are  excitation 
of  electrons  from  the  conduction  or  valence  bands  and  from  the  inner  shells 
of  the  ion  cores  of  the  solid  into  continuum  energy  levels. 
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Four  reports  have  resulted  from  work  performed  under  this  part  of  the 
contract.  These  reports  provide  tabulations  of  inverse  mean  free  paths, 
stopping  power,  ranges  calculated  in  the  continuous-slowing-down  approxima- 
tion (CSDA),  and  straggling  in  ten  solids  for  electrons  of  energies  from  a 
few  eV  through  10  keV.  Also  included  in  these  reports  is  a complete  descrip- 
tion of  the  theoretical  models  and  techniques  employed  in  the  calculations. 

The  solids  studied,  the  report  number,  and  a brief  description  of  the 
material  contained  in  each  report  are  listed  below.  The  references  indicated 
give  complete  information  on  report  title,  number,  and  date  of  publication. 

1.  A1  and  Al^;  AFCRL-TR-75-0583;  Reference  12 

The  interaction  of  electrons  with  the  solids  A1  and  Al^O^  is  described 
based  on  the  electron  gas  model  for  the  conduction  band  electrons  in  Al,  a 
model  insulator  theory  for  the  valence  electrons  in  Al^O^,  and  inner  shell 
ionization  derived  from  atomic,  generalized  oscillator  strengths.  Contribu- 
tions to  the  inverse  mean  free  path  and  stopping  power  from  the  various 
interaction  processes  are  tabulated  for  electron  energies  from  0.5  eV  to 
10  keV  for  Al  and  from  10  eV  to  10  keV  for  Al^O^.  Electron  range  in  the 
continuous-slowing-down  approximation  and  straggling  are  tabulated  for  electron 
energies  from  10  eV  to  10  keV  for  both  materials. 

2.  Si  and  Si09;  RADC-TR-76-1 25;  Reference  13 

The  interaction  of  electrons  with  the  solids  Si  and  SiO^  is  described 
based  on  a "modified  electron  gas  model"  for  the  valence  band  electrons  in 
Si,  a model  semiconductor  theory  for  the  valence  band  electrons  in  S i 0^ » and 
inner  shell  ionization  derived  from  atomic,  generalized  oscillator  strengths. 
Contributions  to  the  inverse  mean  free  path  and  stoppinq  power  from  the 
various  interaction  processes  are  tabulated  for  electron  energies  from 
threshold  4 eV  above  the  Fermi  level)  to  10  keV  for  Si  and  from  10  eV  to 


10  keV  for  SiO^.  Electron  range  in  the  continuous-slowing-down  approximation 
and  straggling  are  tabulated  for  electron  energies  from  10  eV  to  10  keV  for 
both  materials. 

3.  Ni,  Cu,  Ag,  and  Au;  RADC-TR-76-220;  Reference  14 

A statistical  model  is  described  and  employed  to  calculate  inverse  mean 
free  path,  stopping  power,  and  mean  square  energy  loss  for  electrons  of 
energy  from  0.5  eV  to  10  keV  above  the  Fermi  level  in  several  solids.  From 
these  calculations,  electron  range  in  the  continuous-slowing-down  approxima- 
tion and  straggling  are  evaluated  for  electron  energies  from  10  eV  to  10  keV. 
Tables  of  these  quantities  are  presented  for  the  solids  Ni , Cu,  Ag,  and  Au. 
Graphical  displays  of  the  inverse  mean  free  path  and  stopping  power  for 
these  four  solids  are  also  included. 

A brief  description  of  the  use  of  this  statistical  model  in  calculations 
of  electron  inelastic  mean  free  paths  and  energy  loss  is  found  in  Reference  15. 
This  paper  is  also  included  in  this  report  as  Aopendix  H. 

4.  Ge  and  GaAs;  RADC-TR-76-350;  Reference  16 

The  interaction  of  electrons  with  the  solids  Ge  and  GaAs  is  described 
based  on  a model  semiconductor  theory  for  the  valence  band  electrons  and 
inner-shell  ionization  derived  from  classical-binary-col  1 ision  cross  sections. 
Contributions  to  the  inverse  mean  free  path  and  stopping  power  from  the 
various  interaction  processes  are  tabulated  for  electron  energies  from  thres- 
hold (%  2 eV  for  Ge  and  ^ 2.5  eV  for  GaAs)  to  10  keV.  Electron  range  in  the 
continuous-slowing-down  approximation  and  straggling  are  tabulated  for 
electron  energies  from  10  eV  to  10  keV  for  both  materials. 


III.  IMPLICATIONS  FOR  FUTURE  WORK 


The  work  performed  under  this  contract  has  provided  a wealth  of  new  and 
useful  information  concerning  electron  interactions  with  communications 
materials.  However,  the  studies  reported  here  also  raise  many  new  questions 
and  suggest  several  areas  requiring  additional  study. 

One  logical  result  of  our  theoretical  and  experimental  studies  is 

! 

tabulations  of  mean  free  path,  energy  loss,  energy  straggling,  range,  and 
range  straggling  for  electrons  in  materials  of  relevance  to  communications  or 
guidance  systems,  whether  they  be  metals,  semiconductors,  or  insulators.  The 

• models  developed  for  the  tabulations  outlined  in  this  report  are  for  isotropic, 

I 

homogeneous  media.  These  models  may  require  modification  before  they  can  be 
employed  in  the  calculations  necessary  to  produce  the  tabulations  for  organic 
polymers.  In  addition,  the  existing  tabulations  need  to  be  extended  to  lower 
energies.  In  order  to  do  this  it  is  necessary  to  examine,  in  more  detail 
than  is  presently  available,  low-energy  electron  scattering  processes.  Such 
a study  should  consider  electron  elastic  scattering,  which,  though  not 
involved  in  energy  degradation,  must  certainly  affect  transmission  in  a sig- 
nificant manner.  Effects  due  to  the  presence  of  long-  or  short-range 
crystalline  order  in  the  solid,  and  inelastic  (or  quasi-elastic)  scattering 
processes  such  as  phonon  excitation  in  the  "subexcitation"  region  of  energies 
should  be  investigated.  The  effect  of  plasmon  damping  on  low-energy 
electron  mean  free  paths  also  needs  to  be  considered.^ 

The  theoretical  models  involved  in  the  analysis  of  photoelectric  yield 

I 

measurements  to  obtain  electron  attenuation  lengths  require  further  study. 

In  particular  the  shorter  the  values  of  the  attenuation  lengths,  the  more 
important  it  is  to  take  into  account  the  properties  of  the  region  in  the 
vicinity  of  the  surface.  In  the  models  employed  in  the  present  study  planar 
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interfaces  and  homogeneous  media  have  been  assumed.  In  the  UV  and  soft  X-ray 
regions  such  models  have  yielded  electron  attenuation  lengths  comparable  with, 
or  less  than,  the  rms  surface  roughness  height.  Even  assuming  no  other  un- 
certainties, it  is  difficult  to  know  the  significance  of  such  values  in  terms 
of  the  model  assumed  and  the  known  imperfections  of  the  surface  region. 
Finally,  the  theoretical  relationship  between  range,  electron  attenuation 
length,  and  mean  free  path  still  needs  to  be  elucidated. 

Experimentally,  the  problems  encountered  in  this  study  indicate  two 
clear-cut  directions  for  further  investigations:  (a)  to  develop  ways  by 

which  the  extinction  coefficient  of  moderately  to  strongly  absorbing  materials 
might  be  obtained  to  a greater  accuracy  than  is  presently  possible  and 
(b)  to  investigate  ways  of  measuring  electron  attenuation  lengths  which  do 
not  depend  for  their  validity  on  the  homogeniety  of  thin  films  used  in  the 
measurements . 

Specific  projects  which  need  to  be  undertaken  and  which  would  follow 
naturally  from  the  work  performed  under  this  contract  and  its  extensions 
suggested  above  are: 

(i)  Tabulation  of  mean  free  path,  etc.,  for  polystyrene, 

(ii)  Measurement  of  optical  properties  of  other  organic  materials, 
such  as  polyethylene,  and  their  use  in  tabulations  of  mean 
free  path,  etc., 

(iii)  Spatial  distribution  of  Auger  electron-hole  cascades  generated 
by  ionizing  radiations, 

(iv)  Structure  of  heavy  ion  tracks, 
and 

(v)  Slowing-down  spectrum  in  A1  metal. 

Finally,  it  should  be  mentioned  that  several  additional  publications  not 
listed  under  "Publications  and  Reports  Resulting  from  Work  Done  Under  this 


20 


Contract"  (page  22)  will  be  written  based  on  work  done  under  the  contract. 
These  papers  will  report  on  the  slowing-down  spectrum  in  Al^O-j,  the  experimen- 
tal photoelectric  yield  data  as  a function  of  photon  angle  of  incidence,  and 
possibly  give  critiques  on  obtaining  optical  properties  and/or  electron 
attenuation  lengths  from  these  data. 
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ABSTRACT 


We  have  considered  the  generalized  oscillator  strength  (GOS)  corresponding 
to  various  distributions  of  speed  of  a target  electron  in  an  atomic  system  upon 
which  a fast  electron  is  incident.  The  classical  binary  collision  (CBC)  model  is 
employed.  IVe  have  evaluated  the  integral  over  all  energy  transfer  of  the  GOS 
for  each  of  the  distributions  and  find  that  for  large  momentum  transfer  all  of 
those  considered  satisfy  the  Bethe  sum  rule.  At  small  momentum  transfer  all 
fail  to  satisfy  this  rule  due  to  the  fact  that  large  impact  parameter  collisions 
arc  not  described  properly  in  the  CBC  model.  Procedures  for  supplementing 
the  GOS  such  that  experimental  mean  excitation  energies,  polarizabilities,  etc. , 
for  a given  atomic  system  are  reproduced  will  be  discussed  in  a succeeding  paper. 
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I.  Introduction 

Much  progress  has  been  made  in  recent  years  toward  construction  of 
theoretical  cross  sections  for  excitation  and  ionization  of  electrons  from 
stationary  states  in  atomic  systems.  * However,  many  such  approaches  in- 
volve detailed  numerical  solutions  of  the  Schrodinger  or  Dirac  equations  in 
an  assumed  potential  (usually  self-consistent  in  some  sense)  and  require 
separate  studies  for  each  shell  in  an  atom. 

At  a more  empirical  level,  studies  of  the  classical  binary  collision 
(CBC)  model  have  become  popular  following  the  demonstration  by  Gryzinski^ 
that  this  model  has  much  physical  content  and  is  capable  of  giving  quite  rea- 
sonable predictions  of  ionization  cross  sections  for  atomic  and  molecular 
systems  which  may  be  exceedingly  cumbersome  to  treat  on  a more  fundamental 
basis. 

2 3 

The  CBC  model  has  been  applied  * to  the  problem  of  obtaining  estimates 
of  quantities  such  as  stopping  power  for  atomic  systems  for  charged  particles 
with  low  enough  energy  that  the  asymptotic  Bethe-Bloch  formula  is  not  accu- 
rate. These  efforts  suffer  from  the  defect  that  the  CBC  formulas  used  are 
based  on  the  Gryzinski  velocity  distribution  which  is  employed  to  characterize 
the  ground  state  electrons;  a number  of  workers  have  objected  to  this  velocity 
distribution  because  of  its  unphysical,  ad  hoc  character. 

The  authors  have  been  concerned  for  some  time  with  the  study  of  electron 
slowing-down  distributions  in  condensed  matter  at  electron  kinetic  energies 
comparable  with  ionization  energies  in  the  atoms  making  up  the  material.  ^ In 
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order  to  construct  inelastic  differential  cross  sections  for  such  materials  we 
have  studied  the  CBC  model  in  some  detail.  We  have  evolved  certain  simple 
procedures  for  fitting  CBC  model  cross  sections  in  the  large  momentum  trans- 
fer region,  where  they  are  most  reliable,  with  information  on  cross  sections  in 
the  region  of  small  momentum  transfer.  Although  in  essence  semi-empirical, 
our  procedures  result  in  simple  analytic  dif ferential  cross  sections  which  are 
constrained  to  be  accurate  in  the  average  in  certain  respects  described  below. 

They  may  be  employed  to  estimate  charged  particle  inelastic  cross  sections, 
stopping  power,  mean  angular  deflections,  etc. , for  complicated  many-electron 
systems  which  should  be  useful  at  arbitrary  energy. 

Our  approach  is  based  on  sui.-rule  considerations,  but  experimental  data 
on  average  excitation  energies,  polarizabilities , etc.  may  also  be  used  in  this 
connection. 

II.  Classical  Binary  Collision  Model 

Binary-encounter  classical  collision  theory  has  enjoyed  much  attention 

■) 

following  its  application  by  Gryzinski*"  t a number  of  physical  processes  which 

until  then  had  been  treated  mainly  by  quantum- theoretic  methods.  Thomson,  ^ 

7 8 

Thomas  and  Williams  earlier  developer  the  theory  in  some  detail;  however, 
Gryzinski's  work  was  apparently  done  without  knowledge  of  the  efforts  of  these 
authors.  Subsequently  a number  of  mathen  ati  al  < orre<  and  discussions 

9 

of  Gryzinski's  results  took  •!„•  »•  gpev rs  : iut  - n such  theory  is 

correctly  applied  rather  >i  .grec . e aennea;  with  Bom- approximation 


quantum  theoretical  result--  . obtained.  It  is  particularly  valuable  in 

that  it  leads  to  very  convenient  scaling  laws.  Also  it  may  be  applied  to 
systems  which  are  quite  cumbersome  to  treat  using  quantum  methods. 

A.  Differential  C toss  be  tion 

The  basic  doubly-dif  fercrtiai  cro=s  ..  non  on  the  binary  collision 

, . . 6,  7,  8 

model  is  given  by 


p-  = 


(1) 


q 

<£ 


energy  of:  incident  electron 
speed  of  struck  electron  = .2.  £ /m 
momentum  tran-.cr  m collision 
energy  transfer  in  collision 
''ionization  energy"  of  ■ truck  electron  } 


for  spherically  sy  ocities  of  th«  atomic  electron. 

The  law  of  onservation  of  momentum  yields 
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£:,<  £ + ez 


The  quant  r ; ...  tha t thn  cross  section 

r<.  :\  cc  to  • . :•  i L ro:;  ..nd  excited  into 
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B.  Generalized  Oscilla  tor  £ trenirth 

The  first  Born  approximation  in  quantum  mechanics  yields  the  following 

tH 

Expression  for  the  differential  cross  section  for  excitation  of  the  n excited 

\ 

state  in  a system  of  electrons  bombarded  by  a fast  electron  possessing  incident 
energy  E^ ; 


(3; 


for  axially  symmetric  systems.  This  expression  is  valid  in  the  non- relativistic 
limit  and  does  not  include  exchange  effects.  Here 


<:-n|  £e*p(i?.r;)  |o>  (4) 

is  the  matrix  element  of  the  density  operator  evaluated  for  transitions  between 
the  ground  state  with  eigenvector  jo)  and  an  excited  state  with  eigenvector 
}n>»  corresponding  to  the  energy  difference  g - E - E^. 

The  well-known  general  sum  rule*  * may  be  written 

S£lFn('^£)|2  = ^1  Vm  <5> 


where  is  the  number  of  ele  irons  in  t;  system  and  wht  . indicates  a 
sum  over  discrete  final  states  ana  an  inu  -.pal  ov-  • r continuum  final  states. 
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The  differential  1 ciihitoi  ,i.\.  . . ur  transitions  into  the  continuum  may 

be  defined  as 


Ai 

is 

(6) 

In  terms  of  this  quanti  ty 

iTreV*  S-f 

icjiE 

E,£c^ 

(7) 

is  the  doubly  differential  -ross  section  for  excitation  into  continuum  final 
states.  Vriens  studied  the  differential  GOS  for  ionization 

into  continuum  final  states. 

C.  Cliss.ic.al  i'ilo  'vr.t..:  ross  > _uor.s  fo  r W.  rio  • : .-peed 

Distributions  of  the  St  m.  Elciivn 

1.  The  Viria.  Pis ■: ml:  dor.  I.,  this  cast  the  exi  . tetion  energy  is 
■j 

taken  equal  to  £mv,  1 - nt  rah/  . oscillator  stre  h corresponding 

to  the  resulting  dif  ferential  cross  section  may  be  written 


M ^ E,  cf£  AV 
&£  2ire4  <Sd|is 


"Yn  E 


from  Flqs.  (7)  and  (1).  i.  . . (8)  Z , has  U .n  set  oqua*  to  one 


C'.onservatu 


>t  momer 


1 


tend*  \ 
Ktn.r.  I 


a.  u.a 


(9) 
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where  it  is  sufficient  to  consider  E^  >>  (E^,£)  in  evaluating  the  generalized 
oscillator  strength  on  this  model.  If  we  compute  the  integral  of  the  GOS  from 
Eq.  (8)  subject  to  the  restriction  given  in  Eq.  (9)  we  find  that 


1 , f?  > i+  J!  ) 

o , (10, 


where  = q/V^mE^.  In  Fig.  (1)  fv(r))  is  plotted  as  a function  of  tj.  For 
large  enough  momentum  transfer  the  oscillator  strength  defined  according  to 
quantum  mechanics  in  the  Bom  approximation  sums  to  exactly  unity,  as  it 
should.  For  smaller  q the  summed  oscillator  strength  falls  below  this  value 
and  becomes  identically  zero  for  q small.  Thus  one  expects  the  virial  distri- 
bution on  the  classical  binary  collision  model  to  fail  badly  for  small  momentum 

4 9 

transfers  according  to  these  considerations.  Others  ’ have  noted  that  the 
classical  binary  collision  model  is  inadequate  there  because  resonant  energy 
transfers  in  distant  collisions  are  not  included  at  all  in  this  model. 

It  may  be  noted  that  if  Eq.  (1)  is  integrated  over  the  values  of  q allowed 
in  Eq.  (2)  a cross  section  differential  in  energy  is  obtained: 
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1 * £ 


This  integral  can  be  done  analytically.  Setting  £ = ^mv*g,  q = mv0rt>  one 

2 2 

finds  f _ = g ~ g_  when  n < e-  and  f = g + g when  q > e-  Here 
G 1 l Gib  l 

<&.  * -K'+t)  ' I'  - C'+«-)e‘a] 


-it"  6+f)e-‘] 


where  a = 2-q/|  e^r^].  It  may  be  shown  that  f^  1 as  r)-  ® and  - Tpjr 

as  r|  - 0.  Figure  (1)  shows  a plot  of  f vs.  r]  from  the  expression  given 

above  for  the  case  = 1.  It  is  seen  that  the  distribution  chosen  by  Gryzinski 

has  the  desirable  effect  of  yielding  a non-zero  oscillator  strength  as  n -»  0. 

However,  it  is  clear  that  the  value  2/3e^  attained  in  this  limit  has  little 

4 

physical  significance.  Further,  as  Burgess  and  Percival  and  others  have 
noted,  the  Gryzinski  speed  distribution  has  little  basis  in  physical  theory  and 
should  be  regarded  as  empirical  in  character. 

13 

3.  Hydrogenic  Closed  Shell  Speed  Distribution  - Fock  has  shown 
that  the  exact  speed  distribution  corresponding  to  that  of  a closed  shell  of 
any  principal  quantum  number  in  a hydrogenic  atom  is  given  by 


-CfrO  - U 


10 


where  mv2/2  is  the  average  kinetic  energy  o£  electrons  in  the  shell.  This 

O O o 

distribution  is  normalized  such  that  [f(v2)dv2=  !. 

O 

Averaging  Eq.  (1)  over  this  distribution  gives 


o o 


jV 


dV 


^ f W 4le 


(18) 


where  v 


2 m 


. JT , IzL2 

^ m ?<• 


originates  from  the  energy- momentum  con- 


servation conditions  (Eq.  (2))  and  g=  q/®.  The  integration  yields 


i£  __  a'  a/(M* 

Syit  3 f 'c.  / *tj2  j'1 


(19) 


where  the  variables 

\2m £ 

I 

e = e/E 

o 

€•  = € i/E 

<-  o 

p = E./E 
' 1 o 


Ry 


e2/2a 


13.  6 eV 


a = -h2/me2  = 0.529  A 
o 


are  employed. 

The  differential  oscillator  strength  corresponding  to  this  cross  section 
may  be  written 


34 


II 


Kch  . *! 

IT  ’ 3T  (20> 


in  terms  of  the  new  variables.  The  integral  of  this  quantity  may  b<_  done 
analytically;  the  result  is 


where  D = ( ^ - n ) + 4rj  . Figure  1 shows  a plot  of  fc4  as  a function  of 
the  momentum  transfer  variable  choosing  = 1.  It  is  seen  that  the 
velocity  averaging  process  is  successful  in  filling  the  fpe  plane  with  dif- 
ferential oscillator  strength  such  that  f ;h  is  non-zero  everywhere  save  on 
the  n=0  axis,  f CH  tends  asymptotically  to  unity  and  is  slightly  greater  than 

zero  (~  1.  002)  in  the  region  2^  ri  ^ 5-  However,  on  the  rpO  line  there  is 

64  3 4 

difficulty,  since  fCH~  ^ r|  /e-  for  n <<  1.  Bom  approximation  results  with 

various  atomic  models  leave  little  doubt  that  f . should  be  non-zero  there  for 

i 

optically  allowed  transitions. 


III.  Sum  Rule-Constrained  CBC  Cross  Sections 

We  wish  to  retain  the  analytical  simplicity  of  the  CBC  cross  sections, 
but  we  should  prefer  to  use  cross  sections  for  which  the  GOS  satisfies  the 
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fundamental  Bethe  sum  rule1  1 and  possibly  higher-order  sum  rules.  ^ The 
following  simple  procedure  suggests  itself  in  this  connection.  We  note  that 

optical  oscillator  strengths  are  known  experimentally  and  theoretically  for  a 

number  of  atomic  systems  over  a rather  wide  range  of  energies.  ^ Optical 

oscillator  strengths  f (0)  (or  df  (0)/ds  in  the  continuum)  of  an  atom  are  equal 

n 

to  the  GOS  of  that  atom  in  the  limit  of  zero  momentum  transfer;  ie. 


••  't  i^r } 


Suppose  that  f^(0)  and  f(0)  are  known  for  a given  atom  and  that  one  may 
use  CBC  GOS  values  in  the  region  of  Large  r)  for  a given  shell  of  that  atom. 
We  suggest  that  the  following  approximation  to  the  true  GOS  of  a given  shell1 
should  be  useful; 


I C=0  = *«„  W] 


in  the  region  of  discrete  transitions  and 


r,  _ .p  (d)i  ^ 

l <■*  -i  e c 


in  the  continuum.  That  is,  we  propose  an  interpolation  scheme  for  the  GOS 
which  reduces  to  the  known  optical  oscillator  strength  for  small  momentum 
transfers  and  which  goes  over  into  the  CBC  cross  section,  using  a closed  shell 
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hydrogenic  speed  distribution,  tor  Large  momentum  transfer.  There  are  at 

least  two  advantages  to  such  a procedure:  (1)  the  Be  the  sum  rule  for  f and 

— is  satisfied  for  all  q;  (2)  the  cross  section  resulting  from  this  GOS  yields 
de 

a charged  particle  stopping  power  which  reduces  to  the  correct  Bethe  form  for 
large  (but  non- relativistic)  velocities  of  the  charged  particle.  Clearly,  one  may 
choose  other  differential  GOS  functions  than"^  for  use  in  Eq.  (23);  however, 
they  are  plotted  since  the  analytical  simplicity  of  Eq.  (20)  is  quite  appealing. 

Also  this  form  should  be  fairly  accurate  for  inner  shells,  particularly  if  one 
employs  experimental  values  of  the  ionization  energy  E^  of  a given  atomic  shell 

ard  if  one  chooses  the  mean  kinetic  energy  of  electrons  in  that  shell  from  Slater's 

. 17 

rules. 

As  an  indication  of  the  usefulness  of  the  procedure  described  above,  we  have 

plotted  in  Fig.  2 the  exact  dif  fcrential  GOS  of  the  H atom,  df^/de,  the  CBC  GOS 

~ df 

using  Eq.  (20),  and  df/de  from  Eq.  (24)  in  which— (0)  is  taken  to  be  the  exact 

optical  oscillator  strength  of  the  H atom.  All  energies  are  expressed  in  terms 

of  Ry  13.  6 eV  and  all  momenta  are  expressed  in  units  of  (a0)  = (0.  529a.) 

2 2 

The  various  GOS' s are  plotted  as  a function  of  rf  = (a0q)  and  for  several  values 
of  e - fc/Ry.  One  sees  from  this  plot  that  df/de  is  superior  to  dfCi/de  and  is  a 
reasonable  approximation  to  the  exact  GOS  over  the  whole  e-77  plane. 


Summary 

We  have  considered  the  GOS  corresponding  to  various  speed  distributions  in 
the  CBC  model.  We  find  that  the  Gryzin:-.  i distribution,  though  unphysical, 
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yields  a CBC  cross  section  which  satisfies  the  Bethe  sum  rule  better  overall 

than  either  that  derived  from  the  virial  distribution  or  that  derived  from 

the  hydrogenic  closed-shell  distribution.  This  fact  may  account  in  part  for 

the  generally  good  agreement  between  the  Gryzinski  theory-  and  various  ex- 
(4) 

penments. 

In  a succeeding  paper  we  will  present  applications  of  the  approximation 
indicated  in  Equations  (23)  and  (24)  to  the  construction  of  electron  cross  sec- 
tions which  satisfy  the  Bethe  sum  rule  and  yield  mean  excitation  energies 
which  are  constrained  to  agree  with  experiment.  Such  procedures  are  of 
value  in  obtaining  semi--empirical  cross  sections  for  complicated  many- 
electron  systems  to  use  in  the  study  of  electron  slowing-down  in  condensed 
matter. 
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cross  sections  at  small  momentum  transfer,  one  obtains  quite  good  agreement 
with  experimental  ionization  cross  sections  of  the  most  weakly  bound  electrons 
in  O,  N,  Ar,  and  Cl  . Note,  however,  that  he  employs  the  virial  speed  distribu- 
tion in  his  comparisons. 
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FIGURE  CAPTIONS 

Fig.  1.  Integrated  generalized  oscillator  strengths  computed  from  the 
classical  binary  collision  model  plotted  as  a function  of  dimensionless  momentum 

* 

transfer  t> 

Fig.  2.  Generalized  oscillator  strength  for  two  different  values  of  energy 
transfer  c plotted  as  a function  of  2*i(r|  )»  where  n is  momentum  transfer.  The 
quantity  df^/de  is  the  exact  GOS  for  the  H atom,  dfCH/de  is  the  GOS  computed 
from  the  CBC  model  by  averaging  over  the  momentum  distribution  of  the  H atom, 
and  df/de  is  computed  from  Equation  24  of  the  text. 
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APPENDIX  B-l 


Electron  Slowing-  Down  Studies  in  Alo0^ 

Reprinted  from:  Health  Physics  Division  Annual  Progress  Report 

for  Period  Ending  July  31,  1973,  ORNL-4903,  pp.  103-104. 


Because  of  its  importance  as  an  insulator  and  its  ease 
of  fabrication  into  forms  suitable  for  experimental 
study,  aluminum  oxide  was  chosen  for  the  investi- 
gation. The  study  has  been  carried  out,  concomitant 
with  a theoretical  investigation,  using  both  polished  and 
unpolished  polycrystalline  samples.  It  was  necessary  to 
dope  the  AI2Oj  with  a small  amount  of  normal  isotopic 
dysprosium  which  serves  as  the  source  of  primary  beta 
particles  following  activation  by  thermal  neutron  irra- 
diation. 

The  target  disks,  containing  0.5  wt  % dysprosium  as 
the  oxide,  were  fabricated  by  a coprecipitation  tech- 
nique. The  disks  were  ground  to  final  size  using  an 
aluminum  oxide  abrasive.  Sources  obtained  in  this 
manner  exhibited  a density  of  87%  of  the  theoretical 
value.  The  surfaces  of  the  unpolished  disks  were 
covered  with  micropits  and  imperfections  which  had 
dimensions  of  the  order  of  0 001  in.  The  typical  size  of 
the  individual  crystals  appeared  to  be  under  0.005  in. 

Results  of  the  initial  studies  on  the  unpolished 
crystals  were  reported  earlier.  The  electron  flux  inside 

the  material  is  inferred  from  that  measured  in  the 
vacuum  environment  outside  by  dividing  by  a suitable 
barrier  transmission  factor  which  takes  the  form 

AO  = 1 Ve^I(£*£a)  , 


The  earlier  stud.es  were  repeated  using  polished  disks 
to  eliminate  possible  effects  of  surface  roughness  on  the 
experimentally  determined  flux.  A diamond  polishing 
compound  was  used  to  produce  surfaces  with  a flatness 
of  approximately  one  wavelength.  The  results  for  the 
electron  flux  inside  the  specimen,  shown  in  Tig.  10.1, 
exhibit  essentially  the  same  shape  as  the  results  for  the 
unpolished  sources  but  are  everywhere  lower  in  abso- 
lute value.  This  difference  is  thought  to  be  associated 
with  the  uncertainty  in  the  thermal  neutron  flux  in  the 
irradiation  facility  rather  than  in  the  electron  flux  itself. 
The  characteristics  of  the  measured  spectrum  are  similar 
to  those  found  earlier  for  both  metals  and  semicon- 
ductors. The  most  noteworthy  difference  is  the  some- 
what more  rapid  buildup  of  secondary  electrons  at  low 
energies  in  the  case  of  the  aluminum  oxide. 


ORN. -OWC  7J  BO 7? 


where  L is  the  electron  energy  outside  and  I:'A  is  the 
electron  affinity,  equal  to  the  difference  between  the 
work  function  and  the  gap  energy.  There  exists  consid- 
erable uncertainty  in  the  value  of  EA  . A literature 
survey  on  aluminum  oxide  revealed  values  extended 
from  0 for  amorphous  AljOj  to  1.8  eV  for  the 
crystalline  form.  The  electron  affinity  for  our  poly- 
crystalline sources  would  presumably  fall  near  the  lower 
end  of  this  range,  since  this  form  more  nearly  resembles 
the  amorphous  material.  In  our  work  we  have  some- 
what arbitrarily  chosen  an  t'A  value  of  0 5 eV. 


, ' ' •• 

® ' ; t **.. 

t . . ...  ! 

10°  101  10*  103  104 

fcl  t CTPON  INFRGY  (*V1 

Fig.  1 0. 1 . Flectron-siowing-down  spectrum  of  aluminum  oxide. 
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Electron  Slowing-  Down  Studies  in  Al^O^ 

Reprinted  irom:  Health  Physics  Division  Annual  Progress  Report 

for  Period  Ending  July  31,  1974,  ORNL-4979,  p.  147. 


Last  year  we  reported21  measurements  of  the  elec- 
tron slowing-down  flux  in  polished  aluminum  oxide  and 
assumed  that  the  differences  in  flux  values  between 
these  results  and  those  obtained  from  our  earlier  studies 
using  unpolished  sources  resulted  from  uncertainties  in 
the  thermal  neutron  flux  at  the  time  of  irradiation 
rather  than  in  the  electron  llux  itself.  We  have  now 
eliminated  this  uncertainty  through  an  analysis  of  a set 
of  flux  monitors  that  were  irradiated  along  with  the  set 
of  polished  source  disks.  An  accurate  value  for  the 
thermal  flux  was  obtained  from  a cobalt-aluminum 
monitor,  while  corrections  for  the  epithermal  flux  were 
obtained  from  simultaneous  measurements  using  silver- 
aluininuin  monitors.  Uncertainty  in  the  l64Dy  cross 
section  ot  (2700  i 200)  X 10  24  cm2  represents  the 
largest  uncertainty  in  the  determination  of  the  absolute 
secondary  electron  flux 

The  electron  slowing-down  flux  in  aluminum  oxide 
on  the  basis  of  this  new  neutron  flux  determination  is 
shown  in  Fig.  22  12  over  an  energy  range  extending 
from  i .5  eV  to  32  keV.  The  llux  is  slightly  higher  than 
that  reported  previouslv  and  i o flee ts  the  more  careful 
determination  ot  the  2200-m/sec  thermal  neutron  flux 
Also  shown  in  the  figure  is  the  spectrum  calculated  by 
using  Monte  Cailo  techniques  and  employing  inelastic 
cross  sections  obtained  from  a sum-rule-constrained 
classical  binary  collision  model  Delails  of  this  calcula- 
tion are  discussed  in  the  General  Radiation  Physics 
Section 
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ELECTRON  ENERGY  (eV) 

Fig.  22.12.  Hectron-slowing-down  spectrum  of  aluminum 
oxide. 


21  I I Arakawa  el  al.(  Health  Rhys  f) tv.  Anna  I'royr  Rep 
July  .?/.  I'I7J,  ORNI.-490.1,  pp  103  4 
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Reprinted  from:  Health  Physics  Division  Annual  Progress  Report 

for  Period  Ending  June  30,  1975,  ORNI.-5046,  pp.  144-148. 

THEORY  01  THE  LINEAR  RESPONSE  FUNCTION 
OF  A MODEL  INSULATOR2 1 

Inipoitant  piogress  has  boon  made  in  understanding 
tlio  properties  of  metals  through  study  of  the  election 
gas  model  Here  one  ignores  the  detailed  crystalline 
properties  of  the  metal  and  concentrates  on  the 
dynamic  pioperties  of  the  conduction  electrons.  From 
tins  model,  far-reaching  conclusions  have  been  drawn 
about  collective  modes,  screening  of  impurity  charges, 
and  stopping  power,  to  name  only  a few  of  the 
important  phenomena  in  metals  which  have  been 
tieated  on  this  basis. 

There  have  been  a few  calculations  of  the  dielectric 
response  function  of  semiconductors,  semimetals,  and 
insulators.  Most  of  the  recent  papers  are  concerned  with 
the  dielectric  response  function,  e,  of  specific  solids  In 
air  effoit  to  develop  a similarly  simplified  model  for 
insulators  which  retains  useful  physical  information 
about  electronic  properties,  we  have  recently  derived  an 
expression  for  the  linear  response  function  of  a ran- 
domly distributed  assembly  of  molecules.  This  model 
would  be  expected  to  describe  best  the  electronic 
properties  of  an  insulating  liquid. 

The  theoietical  calculation  starts  with  the  polariza- 
tion propagator,  nf  r2 , / : r,  , 0),  describing  the  ground- 
state  expectation  value  for  the  process  of  pair  creation 
at  ( rt , 0)  I ol lowed  by  dest  ruction  of  a pair  at  (r3  t)  It 
is  defined  by 

r(rj,t.  r, , 0)  = iT|  £ t(r; . t)  $(r  2,l)  ^ + (r , , 0 ) 

* ii(r, . 0)1  I sRo)  , (1) 

where  I 'T'o ^ is  the  exact  Heisenberg  state  vector  of  the 
system  (we  consider  only  the  case  7’=0°K),  S2  is  the 
volume  of  the  system,  and  i^(r,  /)  is  the  wave  destruc- 
tion operator  for  electrons  at  (r,  t).  The  wave  operators 
mav  be  expressed  in  occupation  number  representation 
as 

(K  t 

£(r.  /)  21  L “nl(r)dn{0 

I n %.  n | 

un  ore 

* X WkU)2k(0  . (2) 

k 

? I KommuM  sponsored  in  put  by  the  Defense  Nuclear 
Ap*my.  ii lufcr  subtask  IA040,  but  does  not  reflect  endorse- 
ment by  tin*  sponsor. 


46 


145 


where  u„|(r)  is  the  nth  single  particle  "orbital"  on  the 
ion  core  at  position  I and  the  a s ate  the  usual  erection 
operators  We  take  electrons  to  occupy  independent 
orbitals  centered  on  each  of  the  ion  cores,  neglect 
overlap  of  orbitals  on  different  ton  cores,  and  will 
assume  that  there  are  /if  occupied  orbitals  on  each 
core.  With  these  assumptions  we  have  i/„|(r)  un(r  I) 
Although  not  considered  here,  we  could  take  a different 
distribution  of  orbitals  on  different  cores  The  excited 
states  are  taken  to  be  OFW  (orthogonalized  plane  wave) 
orbitals;  that  is, 

e'^  ' y y 

t<k(r)  L L Ormlk)  t/„m(r)  , (3) 

T2  ’ f 2 n < n p m 

where  (r|wm>  = Knm(r)  and  (r  |k>  ’/SI1!2  for  brev- 

ity One  may  easily  check  that  U|((r|  is  orthogonal  to  all 
bound-state  orbitals.  I he  binding  energies  correspond- 
ing to  these  orbitals  will  he  taken  to  fit  experiment 
Since  we  do  not  wish  to  assume  that  the  system 
possesses  any  sort  of  crystal  symmetry,  it  is  not 
appropriate  to  use  Bloch  functions  to  represent  electron 
states. 

The  Fourier  representation  of  n(r2,t;  iq.O)  may  be 
defined  as 

fu3rjtl3tte  ' '■  * ^•r') 


rr( r2 , t,  r,  . 0) 


If  we  abbreviate  fc'q  (2)  as 
^(r,  /)  : £ w„|(r)3„|(Z) 


(0  e ‘^'asi(0)^e  Ssl  , 

(7) 

= *<w*'  2^,(0)  - e‘“*‘ah  , 

where  cos  is  the  cigenenergy  of  an  electron  in  the  sth 
orbital  located  on  any  of  the  ions.  We  also  find  it 
convenient  to  make  the  assumption  that  the  positions 
of  the  ion  cores  are  completely  landom.  Then  we  may 
replace  the  sum  over  I hy  an  integral, 


where  N is  the  number  of  ion  cores  in  the  volume  H. 
With  these  assumptions  we  find 

TfE  I |<p  + q.  |a> 

s P 

£ ( p\n)  (nie  ,q 1 -r|s)  |3  e ’fwp  uj)< 

n 

X 53(q.  + q? ) (8) 

for  t >0  For  the  case  t <0  it  is  only  necessary  to  take 
the  complex  cori|ugate  of  the  time-dependent  exponen- 
tial factor  m Fq.  (8)  We  now  compute  the  frequency 
transform  of  fhe  complete  polarization  propagator  as 


»i|  j .q  j .va 


Ure,UJI  77 C|1 


A' 

= 7>61(4i  f q.-'L  ^ ' <p  + qi  's> 


Fq  (4)  leads  to 

•tq,  ,q2(f ) 

= I e ,(<"  r> 

J.l  K.m  v ,p  vv.q 


2 <p'/t)</il«f  “ll'tls>|1 


x + 

\co(,  uij  to  ir\  Op  la,  t w ii) 


X Hfifrj)  u,im(rj ) u;p(rl  ) irlvq(r,  ) 

X <’l'urn«*t|(f)a1,n1(')«Jp(0)  «„,,(())  | IT„>  ((.) 


To  proceed,  we  need  to  approximate,  we  take  I'k,,) 
to  be  the  state  vector  of  the  nouinteractmg  ground 
state  (|0>),  and  we  take  the  slate  creation  and  annihila- 
tion operators  also  to  have  their  noiunteravting  forms 

Thus 


where  q is  a positive  infinitesimal  Furthermore,  in  case 
the  system  tias  translational  symmetry, 

^qi.q5(f)  -*  5 1 («| i + q2 ) nq|(f); 

in  that  case  the  dielectric  response  function  cq  w is 
given  in  terms  of 

Hq  ,to  /iff  c 7Tq(  t ) 
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as 

‘q.u  - 1 MjAq.w  • (10) 

where  Vq  - 4it/q'  is  the  Fourier  transform  of  the 
Coulomb  interaction  energy.  Note  that  atomic  units  are 
being  used  throughout.  Thus  we  have 

4mV  v 

eq,cj  1 + L 1 pqsl 

' s P 


X 


( 

\U)p  - CJ5  CO  IT] 


I 


tOp  ^ CO 


as  the  dielectric  response  function  for  our  model 
insulator,  where 


„ , „ e dp  * q)  • i 

f"  pqs  = jd3  r 


n'i2 


Wj(r) 


_ r e~ ,p ' 1 1 r 

Zfj\l—-7rujri)fd>r;e''*-'i 


X u„(r2 ) ws(rj ) (12) 

is  the  OPW  form  factor. 

It  is  reasonable  to  look  at  the  case  where  only  a single 
ground-state  orbital  need  be  considered.  This  might  be 
appropriate  when  all  other  orbitals  are  tightly  bound. 
For  an  s orbital  take 


dielectric  response  function  of  a model  insulator  and 
has  been  applied  to  compute  the  differential  inelastic 
mean  free  path  of  a swift  electron  interacting  with 
valence-band  electrons  in  A1203.  We  assume  that  the 
dispersion  relation  for  energetic  electrons  in  the  model 
insulator  is  given  by  u/p  = j}p2 , where  0 = 'J2  if  the 
energetic  electron  has  the  bare  electron  mass,  and  we 
consider  only  one  bound  state.  Instead  of  using  the 
conventional  normalization  for  the  excited-state  OPW 
wave  function,  we  require  that  the  imaginary  part  of  the 
dielectric  response  function,  , satisfy  the  usual  sum 
rule 

CO 

f = 2tt2  n . (1) 

Here  n is  the  density  of  ion  cores,  and  we  continue  to 
work  in  atomic  units.  We  find  the  sum-rule-constrained 
f q2u>  t0  8‘ven  by 

eq2,w  = ™r/p<j A , (2) 

where 


r 


LJ ! 

3 \ [a2  +(d-  p)2V 


[a2  +(/7+P)2  ] 3 


32a1 

(cr+p7  )2(4o2+^:)2 


1 

<*2Hq-p)7 


l 

a2  Hq+p)7 


aN  a properly  normalized  one.  Then 

r.  u foot"'  I 1 

/,pt)0  VTT  { (a2  *-"(q  + p)2J2 

16ur4  | 

(a2  + p 2 )2  (4or  + q2)2  j . (13) 

The  resulting  response  function  is  analytically  quite 
manageable:  it  describes  collective  states  and  screening 
of  impurity  charges  and  at  present  is  being  used  to 
determine  electron  inelastic  mean  free  paths  and  energy 
loss  rates  in  insulators. 


APPI ICATION  OF  rilfc  MODEL  INSULATOR 
THEORY  TO  AUO,:  THE  LOW-ENERGY  END 
OF  lilt  ELECTRON  SLOWING-DOWN 
SPECTRUM22 

The  theory  desenbed  in  the  proceeding  paper  yields 
an  analytical  expression  for  the  imaginary  part  of  the 


I024pqa® 

(a2  +/r  )4  (4a2  +q2  )4 

. (3) 

K + + a2)]  [ -i-,  - 

[ 8a  (4a2+<72)4 

(4) 

p=  [(co  -cog)/0]1/2  . 

(5) 

Here  cog  is  the  excitation  energy  of  the  single  bound 
state.  The  procedure  for  using  this  model  dielectric 
function  to  represent  that  of  a real  insulator  will  be  to 
consider  a and  0 as  adjustable  constants  and  match  the 
optical  dielectric  function  measured  for  that  material. 
In  the  “optical  limit."  q -*  0,  Eq.  (2)  becomes 

(21  _5l2tr  ct1  p3 

3 (6) 
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We  have  tried  matching  e^,2^  as  given  by  Eq.  (6)  to  the 
experimental  values  measured  for  Al203  by  Arakawa 
and  Williams23  Figure  17.11  shows  a comparison  of 
these  results,  where  curve  A represents  the  experimental 
results  and  curve  B is  obtained  from  our  model  with  a = 
0 78,  p = V2 , and  u>g  = 0.3309.  We  used  a density  of 
4.05  for  corundum  and  took  15  valence  electrons  per 
A12Oj  unit.  In  effect  we  have  assumed  that  the 
threshold  for  excitation  of  any  of  the  valence  electrons 
is  9 eV.  We  see  that  the  agreement  is  quite  reasonable. 
We  could  obtain  excellent  agreement  by  taking  a 
combination  of  ground-state  orbitals. 

From  the  values  of  a,  P,  and  ug  determined  in  the 
“optical  limit,"  the  real  part  of  the  dielectric  response 
function,  was  computed  from  e^2^  using  the 

Kramers-Kronig  relations  between  the  real  and  imagi- 
nary parts  of  an  analytic  function.  Figure  17.12  shows  a 
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Fig.  17.11.  The  imaginary  part  of  the  dielectric  response 
function  of  the  model  insulator  in  the  “optical  limit"  (curve  B) 
fitted  to  experimental  results  for  AljOj  (curve  A). 


plot  of  the  energy  loss  function,  Irn  ( l/CqiCJ),  as  a 
function  of  energy  transfer  w for  various  values  of 
momentum  transfer  q from  these  calculations,  where  q 
and  to  are  measured  in  atomic  units  One  sees  a 
prominent  plasmon  resonance  peaking  at  to  ~ I (energy 
transfer  ~ 27.2  cV)  when  q < 1 . As  q increases  to  larger 
values  the  loss  function  peaks  at  values  of  to  which 
increase  and  approach  the  quantity  to  = <72/2,  the 
energy-momentum  relation  of  a free  electron,  as  q 
Thus  our  model  displays  a smooth  transition  in  re- 
sponse from  collective  behavior  for  small  q to  single- 
particle behavior  for  large  q. 

The  differential  inverse  mean  free  path,rip/riw,  for  an 
electron  of  velocity  v interacting  with  a medium  of 
dielectric  response  function  Cq>tJ  is  given  by 


t/p  2 /»V+\A*  5to  (jq  , _j 
du  ~ nvv  vV7Tu 


(7) 


This  integral  has  been  evaluated  numerically  as  a 
function  of  v using  the  dielectric  response  function  of  a 
model  insulator  applied  to  A120j.  The  results  may  be 
fitted  reasonably  well  for  v > 4.7  with  a function  of  the 
form 

rip 

— =/(w)ln(v/v,)  + g(co)  , (8) 

where  /( co)  is  a function  which  describes  mainly 
plasmon  creation,  g(  to)  represents  primarily  interactions 


23.  E.  T.  Arakawa  and  M.  W.  Williams,  J.  Phys.  Client  Solids 
29,  735  (1968). 
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Fig.  I 7.1 2.  Tlie  energy  low  function  a*  a function  of  momentum  and  energy  transfer  from  the  model  insulator  theory. 
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may  thus  prove  lo  have  wide  applicability  to  the  study 
ot  dynamic  electron  properties  in  insulating  substances. 


■ .i  t taiye  ntorncntum  transfers,  and  V!  is 
•mpii  ■ dr:t  • incd  constant  I or  values  of  v < 
4 1 ike  function  .A.  is  stored  in  tabular  form; 

I .„>  e ir  i-rp  at.  in  these  tables  is  carried  out  in  a 
straight). twjiJ  wav  tor  intermediate  values  of  v and  w. 

figure  1 7 | t shows  some  preliminary  results  of  a 
calculation  of  election  slowmg-down  spectra  generated 
in  our  model  insulator  b>  electrons  having  energies  of 
IIX)  and  ;<H)  cV  One  sees  a divergence  in  the  flux  0(F) 
as  the  energv  / approaches  the  threshold  for  excitation 
in  the  insulator  At  the  high-energy  end  of  the  spectrum 
there  are  several  peaks  just  below  the  source  energy 
spaced  at  intervals  of  ~ 2S  eV  These  are  plasmon  loss 
pejks  and  represent  elections  which  are  degraded  by 
the  creation  of  one.  two.  etc.,  pi  asm  oils  after  being 
emitted  from  the  source.  These  calculations  are  being 
extended  to  include  inner-shell  excitations. 

These  calculations  show  that  it  is  possible  to  fit 
experimental  data  on  optical  dielectric  functions  of 
insulating  materials  using  the  model  insulator  theory. 
The  dielectric  response  function  determined  in  this  way 
may  then  he  used  to  calculate  mean  free  paths,  stopping 
powers,  ranges,  slowmg-down  spectra,  etc  This  model 
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for  Period  Ending  June  30,  1976,  ORNL-517 1 , pp.  174-178. 


ELECTRON  SLOWING-DOWN  SPECTRA  IN  Al2Os 

Theoretical  calculations  of  electron  slowing-down 
spectra  for  comparison  with  the  Keplertron  data  require 
detailed  estimates  of  differential  inelastic  cross  sections 
for  excitation  of  valence-band  and  inner-shell  electrons 
for  incident  electron  energies  ranging  from  a few 
electron  volts  to  about  10s  eV.  The  work  reported  here 
involves  the  calculation  of  these  quantities  by  (1)  a 
model  insulator  theory4,5  for  the  excitation  of  elec- 
trons from  the  valence  band  and  (2)  atomic  generalized 
oscillator  strengths  (GOS),  calculated  by  Manson4 * 6  and 
McGuire7 8 *  using  the  Ha rtree -Slater  method,  for  the 
excitation  of  electrons  from  the  aluminum  L shell  and 
oxygen  K shell  and  by  using  hydrogenic  wave  func- 
tions* •’  for  the  excitation  of  electrons  from  the 
aluminum  K shell.  These  cross  sections  have  been 
prepared  for  use  in  a numerical  evaluation  of  a 
Boltzmann-type  slowing-down-cascade  equation.  Con- 
tributions from  Auger  electrons  to  the  slowing-down 
flux  have  been  included  in  a straightforward  manner. 


4.  C.  J.  Tung.  K.  H.  Ritchie,  J.  C.  Ashley,  and  V.  E. 
Anderson,  Inelastic  Interactions  of  Swift  Electrons  in  Solids. 
ORNL/TM-5188  (January  1976). 

5 C.  J.  Tungand  R.  H.  Ritchie,  to  be  published. 

6.  S.  T.  Manson,  Phvs.  Rev.  A 6,  1013  (1972). 

7.  E.  J.  McGuire,  Sandia  research  report  No.  SC-RR-70-406 
(unpublished);  E.  J.  McGuire,  Phys.  Rev.  A 3,  267  (1971). 

8.  M.  C.  Walske,  Ph  D.  thesis. Cornell  University,  1951 ; Phys. 

Rev  101,940(1956). 

9 E.  Merzbacher  and  H.  W Lewis,  Handbuch  der  Physik,  vol. 

XXXIV,  Springer-Verlag,  Berlin,  1958. 


Valence-Band  Excitation 

Since  Al2  03  is  a good  insulator  with  a band  gap  of 
about  9 eV,1 0 we  represent  the  response  of  its  valence 
electrons  using  a model  insulator  theory.  In  the 
previous  annual  report1 1 this  model  was  described  in 
some  detail,  and  the  linear  dielectric  response  function 
derived  from  the  model  was  applied  in  a calculation  of 
the  low-energy  end  of  the  slowing-down  spectrum  in 
Alj  Oj  4 ,s 

Given  the  dielectric  response  function  of  the  model 
insulator,  the  differential  inverse  mean  free  path 
(DIMFP)  is  calculated  from  the  equation4 


t{E,ui)=—=—i  / — Im 

dcj  Trnv  J q q 


where  v = \f2£/m  is  the  velocity  of  the  incident 
electron  and  q±  = (\/2m/h)  (%/f  ± \/E  - hw).  We 
assume  here,  also,  that  the  energy -momentum  relation 
of  the  swift  electron  in  the  solid  does  not  differ 
appreciably  from  that  of  a free  electron  in  a vacuum. 
Figure  30.3  shows  a plot  of  the  quantity  E t{E,  <J)  vs  a> 
for  various  values  of  E,  with  some  quantities  expressed 
in  atomic  units  (a.u.).  In  this  figure,  r(£,  cj)  represents 
the  DIMFP  for  loss  of  energy  co  in  the  excitation  of  an 
electron  from  the  valence  band  by  an  electron  with 
energy  E.  One  notes  when  E > 1 and  u>  > 1 that 
Et  ~ gj'7,  corresponding  to  the  differential  cross  sec- 
tion for  energy  transfer  to  a free  electron.  This  valence 
band  at  9 eV  accounts  for  the  response  of  15  of  the  24 
valence  electrons  per  A1303  unit  in  the  solid.  The 
remaining  nine  electrons  are  assumed  to  form  a second 
tight-binding  level  at  29  eV  and  are  also  described  using 
the  model  insulator  theory. 


Inner-Shell  Excitation 

Theoretical  calculations  of  the  cross  sections  for 
excitation  of  electrons  from  the  inner  shells  of  atoms  in 
a solid  will  be  based  on  two  assumptions:  (1)  that  the 
binding  energies  of  the  tightly  bound  inner  shells  are 
not  changed  substantially  from  the  values  appropriate 
for  isolated  atoms  and  (2)  that  widths  of  the  energy 
bands  corresponding  to  inner  shells  are  narrow  enough 
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Fig.  30.J.  Differential  inverse  mean  free  path,  r(t,  uj).  for 
incident  electrons  of  several  energies  interacting  with  the 
valence  band  (9  eV)  of  ADO j. 

to  be  neglected.  With  these  assumptions,  a theoretical 
determination  of  the  cross  sections  for  the  excitation 
and  ionization  of  inner  shells  may  he  obtained  in  the 
first  Born  approximation  utilizing  atomic  GOS  func- 
tions for  such  inner  shells.  The  details  of  our  calcula- 
tions of  DIMFPs  for  inner-shell  excitation  have  been 
presented  earlier  and  need  not  be  repeated  here  1 213 

Auger  Electrons 

Auger  electrons  are  produced  when  vacancies  in  a 
given  shell  are  filled  by  electrons  from  a higher  shell. 


10  F T Arakawa  .ind  M R Williams,  J Phys  them  Solids 
29.  735  (|UM| 

11  R H Ritchie  et  .it  . Health  f'hvs  Div  -1  nnu  Prog  Rep 
June  JO.  It’S  ORNL  5046.pp  144- 48 

12  Ibid  .pp  154  -55. 

13  J C Ashley,  C 1 Tung,  ami  K II  Ritchie.  IEEE  Trans 
Noel  Sci  NS-22.  2533  1 1975) 


The  Augei  transition  lesults  in  a two-electron  process. 
One  election  fills  the  hole  in  the  inner  shell,  while  the 
excess  of  binding  energy  is  transferred  to  a second 
election  causing  Us  emission  as  an  Augei  election  l’u 
oh  I a n i Augei  cascade  I r ansi  lions  In  Alj  Ui , we  make  l hr 
following  mtsuinptlons  II)  tlml  llic  Augei  liuriMtloiia 
on  nr  only  on  u given  ion  (elllirr  aluminum  or  oxygen) 
and  not  between  ions  and  (2)  that  the  Auger  transitions 
take  place  only  from  levels  immediately  above  those  in 
which  holes  exist.14  With  these  assumptions,  the 
iransnion  to  a vacancy  in  the  aluminum  K shell  can  be 
from  either  the  aluminum  L , shell  (2s)  or  the  alumi- 
num Le  , shell  (2p).  We  have  used  the  transition  from 
the  aluminum  Z,2  3 shell  to  a vacancy  in  the  aluminum 
K shell  because  selection  rules  tend  to  favor  this 
transition.  Also,  calculation  of  the  Auger  electron 
kinetic  energy  for  each  transition  and  the  number  of 
Auger  electrons  crealed  following  the  filling  of  a single 
hole  in  a given  shell  is  straightforward.  For  example,  the 
filling  of  a vacancy  in  the  aluminum  K shell  from  an 
electron  in  the  aluminum  2-j.j  shell  will  create  an 
Auger  electron  having  energy  equal  to  hcoA,A 
2tkjAI"f  :2'3  — 1384  eV,  where  heu^  is  the  binding 
energy  of  the  indicated  level,  and  will  p've  rise  to  two 
holes  in  the  valence  band  and  an  Augei  electron  with 
energy  a 63  eV. 

Electron  Slowing-Down  Cascade 
Equation  and  Its  Solution 

Spencer  and  Fano1  5 and  Spencer  and  Attix16  have 
pven  a theoretical  treatment  of  electron  slowing-dowm 
spectra  at  energies  high  enough  that  detailed  electronic 
structure  of  the  material  is  not  important.  In  their 
treatment,  the  M6Uer  cross  section  for  electron-elec- 
tron scattering  was  employed  to  compute  the  delta-ray 
flux  generated  by  electrons  interacting  with  the  medi- 
um. This  delta-rav  flux,  when  added  to  the  flux  of 
electrons  slowing  down  from  higher  energies,  gives  the 
total  flux  in  the  medium  at  arbitrary  energy  The 
Spencer-Fano  approach,  as  oripnally  employed,  is 
strictly  valid  only  for  electron  enerpes  that  are  large 
compared  with  the  binding  energies  of  the  most  tightly 
bound  electrons  in  the  medium.  Calculations  at  enerpes 
comparable  with  transition  energies  from  inner  shells 
require  not  only  the  detailed  cross  sections  but  also  a 
different  approach  to  the  solution  of  the  Boltzmann 
equation.1  7,1  8 

Here  our  interest  is  in  the  energy  dependence  of  the 
flux  of  electrons  generated  in  a solid  irrespective  of 
their  spatial  distribution  The  Boltzmann-tvpe  electron- 


52 


number-balance  integral  equation  lor  the  flux  may  be 
written  easily.  Consider  a vers  laige  medium  containing 
an  election  source  density  that  is  uniform  in  spam  .d 
constant  in  tune  and  produces  .S’tf'l  df  elections  in  the 
interval  of  energy  between  t ano  £'+. IE  per  unit 
volume  The  electrons  released  trom  this  source  lose 
energy  in  collisions  with  othei  electrons  m the  medium 
In  the  steady  state,  the  number  of  elections  that  enti  a 
given  interval  of  energy  between  E and  t + dt 
equal  the  number  scattered  out  Tlius.  ;lie  elect/rn 
slowing-down  cascade  equation  i:.  terms  of  ihe  flux 
<p(t’)  may  be  written  as 

0f£")  f dui  t(E.  oj)  = f du>  <HE '+  oe>  |:|F ' + w. 

J o J o 

+ Tj (£"  + w,£')] +S(i‘)  . (2) 

where  t(E.  to)  is  the  D1MFP  that  an  electron  with 
energy  E may  lose  energy  in  the  interval  between  u and 
co  + dco  and  t s(£,  E)  is  die  D1MFP  that  an  electron 
with  energy  E will  produce  a secondary  electron  in  the 
interval  of  energies  between  E‘  and  E'  r dE'  This 
equation  may  be  solved  by  standard  numerical  tech 
niques,4  given  t and  ts  for  a solid  of  interest 
Numerical  evaluation  of  <ME)  from  tq.  (2)  ovei  die 
energy  range  of  the  Keplertron  data19  would  require 
prohibitively  large  amounts  of  computer  time,  since  it  is 
necessary  to  use  an  energy'  mesh  size  comparable  with 
or  smaller  than  a few  electron  volts  (the  mag.iitude  of 
valence-band  excitation  energy)  in  order  to  solve  Lq. 
(2)  accurately  by  direct  methods  Hence,  we  have 
assumed  that  the  Spencer-Fano  (SF)  flux1'  vs ( (E)  is 
sufficiently  accurate  down  to  Em  ~~  10  keV  We  have 
determined  an  equivalent  source  density  Ss Ft / ) below 
Em  by  computing  an  equivalent  first-collision  density  in 
th'x  range  using  4>sy(E)  and  D1MFI’  functions  trom  die 
model  insulator  theory  and  from  the  inner-shcli  tioss 


14  D Penn.  Rhys  Rev  128,2093(1962) 

15  L V Spencer  and  U Fano.  Rhvs  Ret  93.  1 1 :2 '.  19541 

16  L V Spencer  and  I It  Attix  RaJtat  Ret  3.  239 
H955) 

17  l C hmeuon.k  D Birkholf,  V t Anderson.  and  R It 
Kilchie.  Rhys  Rev  R 7.  t 798  (197  3) 

18  R H Ritchie,  ( I Tunir.  V'  I Andre-.  I l 

Ashley.  Radial  Res  64.  181  0975) 

19  E T Arakau.i  cl  al  . Health  Rh  is  IJi  i Art  m Rr.  >g  Ref 
July  31.  1974  ORNI  -4979.  p 147 


section  calculations  Figure  30.4  shows  a plot  ot  the 
Spencer  Fano  dux  for  energies  above  1 keV  due  to 
beta  rays  front  neutron-activated  16“Dy  uniformly 
distributed  m the  Al20s  sample.  The  flux  has  been 
MirmaJi/ed  to  one  electron  created  per  cubic  centi- 
meter, so  that  the  units  of  this  normalized  llux  are 
entimeters  per  electron  volt.  In  Fig.  30. S,  we  show  the 
. urce  density  Ssf(E)  calculated  using  Eq  (2)  with  the 
Spencer  Fario  dux  <psl  (£)  and  D1MFP  functions  tor 
A1.O3.  At  high  energies,  most  of  the  contribution  to 
the  equivalent  source  density,  Sst-iE),  is  from  degraded 
Sjieneer  Fano  llux  electrons  from  higher  energies.  As 
lectrort  energy  decreases,  the  source  density  increases, 
nee  more  and  more  secondary  electrons  are  generated. 
The  number  ot  secondary  electrons  generated  m the 
medium  increases  very  rapidly  as  electron  energies 
approach  the  threshold  energy.  In  the  bottom  left 
■met  oi  Fig.  30.5,  a separate  plot  of  the  function  for 
r riles  very  close  to  Em  shows  the  general  behavior  of 
the  function  in  that  energy  region. 

Further  details  of  the  calculation  of  the  slowing-down 
spectrum,  given  the  equivalent  source  density,  are 
rather  lengthy  and  will  not  be  repeated  here  (see  ref.  4). 
Figure  30.6  shows  a plot  ot  a measured  electron 
owing-down  flux  in  AlvOj  due  to  beta  tays  from 
neutron-activated  164Dy  uniformly  distributed  in  the 
ample  and  the  result  ol  our  corresponding  calculations 
Isolid  line).  For  comparison,  we  show  a Monte  Carlo- 
computed  spectrum18  using  semiclassical  binary  colli- 
sion cross  sections.  Ihe  theoretical  flux  seems  to 
onlirm  ihe  Monte  Carlo  calculations  but  shows  some 
discrepancy  between  experiment  and  theory  in  the 
energy  range  from  40  to  2000  eV.  Efforts  to  under- 
stand these  discrepancies  are  continuing.  However, 
general  agreement  between  theory  and  experiment  in  a 
broad  energy  range  is  clear 
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Fig.  30.6.  Electron  slowing-down  spectrum  in  Ai203.  The 
solid  curve  shows  the  result  of  the  present  calculation.  The 
points  indicate  experimental  values.  The  histogram  shows  the 
result  of  a theoretical  Monte  Carlo  calculation. 
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Fig.  30.5.  First-collision  electron  source  density . .Vs, [1  6),  in 
AlsOj,  calculated  from  differential  inverse  mean  tree  path 
functions  derived  in  this  work,  and  from  the  Spencer -Fano 
flux  for  electron  energy  above  10  keV. 
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6.13.  Electron  Attenuation  Lengths  in  Carbon  hlms1 * *) 

E.  T Arakawa,  M.  W.  Williams,  R N Hamm  and  R 0.  Birkhoff 

Health  Physics  Division,  Oak  Ridge  National  Latxi'^iory.  Oak  Ridge.  Tennessee  37830  USA 

A J.  Bra  ndmeier 

Department  oi  Physics,  Southern  Illinois  University.  Edwardsville,  Illinois  620?b 

A method  lias  been  developed  Tor  obtaining  electron  attenuation  lengths  in  solids  fiom  measurements  of’ 
photoelectric  yield  as  a 1 unci  ion  of  photon  angle  of  incidence  1 1 j Conventionally,  election  mean  tree  paths, 
or  the  directly  measured  quantities  electron  attenuation  lengths,  are  determined  from  electron  transmission 
measuiements  through  thin  films.  However,  when  electron  attenuation  lengths  are  short,  it  is  difficult  to 
prepare  sufficiently  thin  homogeneous  films  for  these  thin  film  measurements.  Also  when  energetic 
elections  interact  with  a solid,  the  points  of  origin  and  eneigies  ot  secondaty  electrons  are  not  well  defined , 
and  it  is  difficult  to  extiact  a meaningful  electron  attenuation  length  from  such  measurements  The  method 
currently  developed  for  obtaining  electron  attenuation  lengths  avoids  these  difficulties  Thick  films  of  the 
solid  under  investigation  are  used,  avoiding  the  problems  associated  with  making  very  thin  homogeneous 
films  Furthermore,  the  quantity  analyzed  is  the  photoelectric  yield  as  a function  of  the  photon  angle  of 
incidence  lelativc  to  the  photoelectric  yield  for  normal  incidence,  so  that  absolute  measurements  are  not 
required.  Absolute  photoelectric  yields  are  difficult  t<»  nioaone  accurately,  hut  relative  values  are  easy  to 
obtain  since  virtually  all  phutuelectrons  emitted  in  any  direction  from  the  photoemilter  can  he  collected  if 
d sufficiently  high  poteniiul  difference  jx  applied. 

I he  expei  imental  apparatus  has  been  described  previously  A grazing  incidence  scanning  monochtomator  ( 2 ] 
is  optically  coupled  to  a stainless-steel  high  vacuum  chambei  [3).  The  essential  features  of  this  chamber  are 
a sample  holder  at  the  centei  ot  a cylindrical  electron  collector  1 lie  sample  holder  can  he  rotated  through 
360°,  so  that  the  angle  of  incidence  of  the  photons  can  he  varied  1 he  cylindncal  election  collector  allows 
collection  of  virtually  all  photoelectrons  emitted  in  am  direction  Thick  carbon  films  were  obtained  by 
vapor  deposition  from  a carbon  arc,  in  vacuum,  directly  onto  a clean  glass  slide  Alter  insertion  of  a carbon 
film,  the  experimental  chamber  was  pumped  down  to  ~ 2 X 10  7 torr. 
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f or  each  photon  energy  the  total  photoelectric  yield  was  obtained  in  terms  ot  tlie  current  going  from  fhe 
sample  to  t!ie  election  collector,  as  a t unction  of  photon  angle  of  incidence,  0 The  yield  at  angle  0,  Y(h) 
relative  to  the  yield  for  normal  incidence,  Y(0).  was  analy/ed  a>  a function  of  0 as  described  in  Kef  j 1 1 
For  a givei-  photon  energy  the  ratio  Y(P)/Y(0)  (ipO  slu*w  ^ .» single  maximum  when  plotted  versus  0 j-. 
shown  in  Fig  l it  the  t atm  is  intimtely  thick  ami  it  Ms  letiactive  index,  u.  is  less  than  un.'v  hoi  a small 
extinction  coefficient,  k.  ( 0 l)  these  cuives » an  lv  uuly/ed  to  yield  it  and  eithei  k oi  1 wheic  I i>  tie* 
electron  attenuation  length  tor  the  photoexuted  elections  I he  position  ol  the  maximum  of  G(0)  vs  0 is  m 
the  vicinity  of  0C  where  (\  sin  1 n and  hence  yields  the  value  of  n while  the  magnitude  ot  G(flc)  i* 
sensitive  to  k and  I In  the  present  case  we  obtained  values  ot  k independently  from  measurements  of 
reflectance  as  a function  of  photon  angle  of  incidence  and  then  used  a two  parameter  least  squares  fitting 
procedure  to  obtain  n and  L from  (i(0)  vs  0.  f ig  1 shows  representative  data  points  (shown  by  ) for  a 
photon  energy  ot  51  87  eV.  Using  an  independently  measured  value  of  k 0 04  ,i  least  squares  fit  of  this 
data  to  the  theory  (shown  by  •)  gave  n 0.89  and  L - 2.4  A.  Since  (he  reflectance  vs.  photon  angle  of 
incidence  measurements  yield  both  n and  k values,  we  can  make  a consistency  check  bv  comparing  the 
values  of  n obtained  from  the  reflectance  and  from  the  photoelectric  yield  measurements 
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Fig.  2 

Electron  attenuation  length  loi 
photoexuted  electrons  in  carbon  \s 
incident  photon  energy 


The  values  obtained  for  the  electron  attenuation  length  of  photoexcited  electrons  in  carbon  films  are  shown 
in  Fig  2 for  photon  energies  from  20  to  04  eV.  This  is  the  range  of  energies  for  which  the  conditions  n ^ 1 
and  k <•  0 1 are  satisfied  for  carbon. 

Note  added  m proof  Mote  recent  experiments  on  the  total  photoelectric  yield  as  a function  of  carbon  film 
thickness  lead  us  to  suspect  that  the  values  of  the  electron  attenuation  lengths  at  the  lower  photon  energies 
( 20  30  eV)  should  be  considerably  less  than  shown  in  Fig.  2 
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ELECTRON  ATTENUATION  LENGTHS 
IN  CARBON  FILMS25 

Electron  attenuation  lengths.  A,  of  photoexcited 
electrons  in  thick  vapor  deposited  carbon  films  have 
been  reported  previously2  6 over  the  range  of  photon 
energies  from  20  to  64  eV.  These  values  were  derived 
from  measurements  of  the  total  photoelectric  yield. 
Y(0 ),  as  a function  of  photon  angle  of  incidence.  0 For 
a small  extinction  coefficient,  k (<0.1 ),  and  a refractive 
index,  m less  than  unity,  the  ratio  Y(8)IY( 0)  - (i'(O)vs 
0 curves  were  analyzed  to  yield  n and  A.  For  each 
photon  energy  the  position  of  the  maximum  of  Off*)  vs 
6,  which  occurs  in  the  vicinity  of  0 , where  0 = sm  1 n, 
yields  a value  of  n while  the  magnitude  of  G (0 .)  is 
sensitive  to  boih  k and  A We  measured  k independently 
and  hence  obtained  values  of  ii  and  A by  this  method. 
As  a consistency  check,  we  compared  the  values  of  n 
obtained  from  reflectance  and  from  photoelectric  yield 
measurements  on  similar  samples  and  obtained  very 
good  agreement.  However,  even  though  this  consistency 
check  was  satisfactory,  we  still  questioned  the  high 
values  of  / ( — 1 AO  to  50  A)  obtained  in  the  photon 
energy  region  of  20  to  30  eV  We  thus  decided  to 
obtain  values  of  A from  measurements  of  the  photo- 
electric yield,  ) (r).  as  a function  ot  carbon  film 
thickness,  I.  A wedge-shaped  carbon  film  was  deposited 
on  a glass  slide,  and  the  metrical  thickness  as  a 1 unction 
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of  position  on  the  slide  was  obtained  fiotn  a combina- 
tion of  tiansmittance  and  independent  interferometric 
measurements.  The  experimental  apparatus  described 
previously'*  was  used.  In  the  same  way  as  for  the  ) (0) 
vs  0 measurements,  the  carbon  films  weie  deposited  in 
vacuum  and  exposed  to  the  atmosphere  on  being 
transferred  to  die  apparatus  used  tor  photoelectric 
measurements;  then  the  values  of  Y(t)  were  obtained 
with  the  experimental  chamber  pumped  dow  n to  ~ 2 X 
lCT7  torr.  For  each  constant  incident  photon  energy,  a 
Y(t)  vs  t curve  was  plotted.  Some  representative  curves 
of  Y(t)  vs  t for  carbon  on  glass  are  shown  in  Fig.  1 6.8. 
It  was  found  that  for  each  photon  energy  from  10  to  64 
eV,  there  was  a rapid  increase  in  F(r)  in  the  vicinity  of  1 
= (10  to  15)  A with  saturation  of  Y(t)  above  this 
thickness.  The  complete  analysis  of  photoelectric  yield 

ORNL-DV.'G  75  6468 


Fig.  16. ft.  Photoelectric  yield  vs  carbon  film  thickness  for 
constant  incident  photon  energies. 


measurements  from  such  thin  films  is  very  complicated, 
requiring  a knowledge  of  the  homogeneity  of  the  films 
and  their  effective  optical  piopertics  in  addition  to  a 
knowledge  of  the  factois  affecting  the  electron  trans- 
port properties.  Although  a complete  analysis  of  our 
1 (i)  vs  t data  for  carbon  lias  not  yet  been  attempted, 
the  values  of  / for  pliotoexcited  electrons  in  carbon 
appear  to  be  much  smaller  in  the  region  of  20  to  40  eV 
than  the  ones  obtained2"  from  the  observations  of  ) (0 ) 
vs  0.  A theoretical  analysis  of  the  )'(0)  vs  0 method  has 
since  been  undertaken  which  showed  that  although 
reliable  values  of  n can  be  obtained  from  <7(0)  vs  0.  to 
obtain  /.  to  il  A requires  that  k be  known  to  within 
~1T>.  The  uncertainty  in  the  independently  measured  k 
values  was  much  greater  than  this,  loading  to  very  large 
uncertainties  in  the  calculated  /.  values. 

Electron  attenuation  lengths  arc  known  to  be  small  m 
the  region  of  10  to  100  eV,  and  hence  then  accurate 
determination  is  difficult  At  the  same  time,  the  values 
of  these  electron  attenuation  lengths  in  both  metals  and 
insulators  are  of  gieat  interest.  We  thus  regard  this  as  a 
continuing  project.  For  the  F(0)  vs  0 method,  we  pljn 
to  investigate  vaiious  techniques  whereby  we  can  obtain 
extinction  coefficients,  k,  to  a higher  accuracy  than 
presently  attainable  since  we  have  shown  that  the 
uncertainties  associated  with  the  values  obtained  for  /. 
.ne  directly  iclated  to  the  uncertainties  in  I lie  input 
values  of  k For  the  ) (r ) vs  t method,  we  plan  to 
prepaie  and  measure  films  of  gold  in  ultrahigh  vacuum 
for  comparison  with  preliminary  data  already  obtained 
on  gold  wedges  piepaicd  in  vacuum  but  then  exposed 
to  air  on  transfer  to  the  apparatus  used  to  measure  >'(/) 
vs  t.  Gold  is  simpler  to  work  with  than  carbon  and 
should  provide  a check  on  the  reliability  of  the 
experimental  >'(<)  vs  I curves  obtained  for  carbon. 
Furthermore,  experimental  and  theoretical  studies  of 
the  electron  attenuation  lengths  and  the  related  elec- 
tron mean  free  paths  for  gold  have  been  reported  in  the 
liteiature.  Data  on  the  effective  optical  properties  of 
gold  as  a function  of  film  thickness  are  also  available  for 
evaporated  gold  films  ptepated  under  various  condi 
tions.  A complete  analysis  of  the  Y(t)  vs  t data  for  gold 
should  thus  be  possible,  providing  a check  on  the 
reliability  of  this  method.  Ultimately,  we  hope  to  get 
completely  consistent  and  reliable  values  of  /.  in  both 
metals  and  insulators. 
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Opr  ual  tiuKt  tom;  of  po  1 vstyrene  in  th  • torn  of  thin  f i I ms  were 
determined  lor  photon  energies  between  ().(>  and  82  <.■  V l ro.u  t rattsni ion 
measurements . The  results  lot  I:,  the  extinction  coeffic  ient , were* 
combined  with  previous  experimental  results  in  the  soit  and  hard  X-ray 
regions  up  to  8050  eV.  Analyses  were  made  on  several  sura  rules  for  the 
optical  functions  in  this  unusally  wide  energy  range,  including  a sum 
rule  for  the  refractive  index,  n,  derived  recently  by  Altarelli  et  al. 
Redistribution  of  the  oscillator  strength  corresponding  to  2.17,  of  the 
total  electrons  was  found  between  the  valence  and  core  excitations  of  carbon. 
Using  the  complex  dielectric  function  and  the  energy  loss  function  obtained, 
the  average  photoexcitation  energy  and  the  average  energy  loss  for  fast 
charged  particles  over  the  entire  oscillator  strength  distribution  were 
evaluated  to  be  25.1  and  36.8  eV,  respectively. 

INTRODUCTION 

The  primary  purpose  of  this  paper  is  to  examine  in  detail  the 
response  of  electrons  in  a typical  solid  hydrocarbon  to  electromagnetic 
radiation  over  essentially  the  whole  frequency  range.  This  is  achieved 
through  sum  rule  calculations  which  are  employed  to  determine  the  way 
in  which  the  oscillator  strength  converges.  In  the  past  such  calculation-; 
and  their  interpretation  have  suffered  to  a greater  or  lesser  extent 
because  ihe  experimental  data  were  generally  unavailable  over  a wide 
enough  energy  range  to  give  the  desired  certainty  in  the  required 
extrapolations  to  infinite  energy.  To  our  knowledge  there  are  just 
two  previous  studies  performed  on  experimental  data  which  extends 
over  a sufficiently  large  energy  range.  V.’u  know  of  no  such  analysis 
for  any  hydrocarbon. 
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In  this  study,  we  report  the  optical  functions  of  polystyrene 
-(CgH^)-,  in  the  form  of  thin  solid  films  over  the  range  of  photon 
energies  from  0.6  to  82  eV.  The  results  for  k,  the  extinction  coeffi- 
cient, are  combined  with  the  previous  experimental  data  from  30  to  1550 

2 3 

eV  taken  by  Lukirskii  et  al.  and  those  from  1490  to  8050  eV  by  Nordfors. 

The  combined  spectrum  of  k covers  substantially  all  the  electronic 
excitations  (more  than  99.9%  of  the  total  oscillator  strength),  includ- 
ing the  carbon  core  excitations  starting  at  around  282  eV.  Despite  the 
recent  progress  in  synchrotron  radiation  techniques,  which  make  it 
possible  to  obtain  optical  spectra  from  the  infrared  to  the  hard  X-ray 
region  using  a continuum  photon  source,  similar  data  for  any  kind  of 
material  are  still  quite  rare. 

4 

In  addition  to  the  well-known  f-sum  rules,  a detailed  analysis  is 

performed  in  this  study  on  the  refractive  index,  n,  at  high  energies  in 

order  to  employ  a new  sum  rule  for  n which  has  been  derived  recently  by 
5 

Altarelli  et  al.  using  the  superconvergence  theorem  in  high-energy 
physics.  Since,  as  will  be  shown,  the  convergence  of  this  sum  rule  is 
quite  slow,  it  can  he  demonstrated  and/or  used  only  when  the  n values 
are  available  up  to  the  extremely  high  energy  region. 


EXPERIMENTAL 

Optical  functions  of  polystyrene  films  between  0.6  and  82  eV  were 
determined  by  the  transmission  method  , details  of  which  have  been 
described  previously.^  In  this  method,  the  quantities  determined 
experimentally  are  the  transmittance,  T,  of  a free-standing  film  as  a 
function  of  photon  energy  over  the  whole  range  of  photon  energies,  and 
the  film  thickness,  d.  In  addition  the  values  of  n are  measured 
directly  over  a limited  energy  region  where  the  films  are  transparent. 


To  derive  the  optical  functions  n and  k at  each  energy,  the 


3 


4 


Kraners-Kron Lg  relation  between  n anil  k is  then  utilized.  An  iterative 
computational  procedure  is  carried  out  between  the  Kramers-Kronig 
relation  and  the  explicit  expression  for  transmittance  which  includes 
the  correction  factor  for  reflection  effects  at  the  film  boundaries, 
until  the  values  of  n and  k converge,  n and  k are  thus  determined  self- 
consistently  so  that  they  reproduce  the  observed  value  of  transmittance. 

Self-supporting  films  of  polystyrene  were  prepared  by  the  same 
method  as  described  previously.^  Uniformity  of  film  thickness  was  checked 
by  measuring  the  transmittance  at  a suitable  wavelength  while  scanning  the 
whole  film  area.  Determination  of  the  film  thickness,  d,  was  made  by 
observation  of  transmission  maxima  and  minima  due  to  interference  in  the 
transparent  region.  For  the  thinnest  films,  interference  patterns  could 

not  be  used  as  even  the  first  minimum  did  not  occur  in  the  non-aboscbing 

o 

region  below  4.4  eV  (above  ^2800A) . In  such  cases  the  thickness  was 
determined  directly  from  the  explicit  expression  for  the  transmittance 
in  the  transparent  (k=0)  region,*^ 
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Ki  tt.iv  t ivi  indues,  11,  between  2.0  and  4.4  eV,  where*  polystyrene  is 
transparent,  were  determined  by  the  critical-angle,  Brewster  angle,  and 
interference  methods  using  separately  prepared  films.  All  these  methods 
gave  consistent  results  for  n to  within  an  error  of  +0.5%.  The  average 
values  were  used  in  the  thickness  determination. 

g 

As  is  well-known,  polystyrene  shows  a weak  absorption  around  4.8  eV 

O 

(a-2600  A),  which  is  a remnant  of  the  forbidden  absorption  in  benzene  at 
about  the  same  energy.  Since  this  absorption  was  too  weak  to  be  observed 
even  in  the  thickest  of  the  films  used  in  the  observations  described  so  far, 
an  additional  transmission  measurement  was  performed  in  the  vicinity  of  this 
band  on  a 4.50p-thick  film,  and  the  extinction  coefficient  k for  this  band 
determined.  The  measurments  were  extended  on  this  thick  film  down  to  0.6  eV , 
where  no  further  absorption  was  observed.  This  thick  film  was  also  used  to 

determine  the  density  of  the  polystyrene  films.  A knoi/n  area  of  the  film 

3 

was  cut  out  and  weighed.  The  density  was  found  to  be  (1.05±0.02)  g/cra  . 

g 

This  value  agrees  with  the  literature  values  for  bulk  polystyrene:  (l .042-+ 

3 

1.065)g/cm  . The  tacticity  of  polystyrene  in  our  films  was  not  checked. 

It  does  not,  however,  seem  to  affect  the  results  in  the  present  spectral  region. 
Errors  in  the  results  obtained  for  (n-1)  and  k were  estimated  to  be  less  th3n 
+2%,  except  for  the  region  above  60  eV. 

The*  measurements  described  above  were  made  with  a Cary  recit'd  * ; 
spec t rophotom-.'ter  (model  14)  in  the  region  0.6-2  eV,  a Seya-Namioki 
monochromator  (McPherson  model  215)  in  the  region  2+15  eV,  and  a grazing- 
incidence  monochromator  (McPherson  model  247)  in  the  region  11-S2  eV . 

Above  7 .5  eV  the  measurements  were  made  with  I ine  sources  (hydrogen  rttso.il  ir'e 
anJ  a condensed  air  spark)  while  below  7.5  eV  continuum  sources  wen*  employed. 
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RESULTS  AND  DISCUSSION 

1.  The  extinction  coefficient  k. 

The  results  for  k between  0.6  and  82  eV  are  presented  in  Fig.  1 on 
a logarithmic  plot.  The  spectrum  exhibits  structure  below  10  eV  due  to 
molecular  electronic  excitations.  Above  10  eV,  the  spectrum  is  essentially 
structureless,  except  for  a broad  absorption  peak  at  around  14  eV.  The 
absorption  decreases  rapidly  and  monotonically  above  this  peak,  and  can  be 
described  by  a straight  line  above  20  eV  on  a logarithmic  plot.  The 
best  fit  of  the  present  k values  between  20  and  60  eV  to  a straight  line 
gave  the  expression 


k(E)  = 7.81  x 103  E3’23,  (2) 

where  E is  the  photon  energy  in  eV.  Deviation  found  above  60  eV  of  the 
k values  from  the  straight  line  fit  seems  to  be  due  to  experimental  errors 

in  the  transmission  measurements  in  that  region,  where  the  transmit tances 
of  thin  films  were  quite  high. 

Also  shown  in  Fig.  1 are  the  combined  results  for  k obtained  by 
2 3 

Lukirskii  et  al.  and  by  Nordfors  extending  to  the  hard  X-ray  region. 

o 

These  results  wire  originally  given  by  the  expressions'1 


2 T 1 

p(\)  = 0.40.1 


(3) 


o o 

for  44  A < X < 410  A and 


u(A> 


2.4012 


.68 


(4) 
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for  8 A •-  ' *>44  A,  wfu-re  p is  Che  linear  absorption  coefficient  in 


cm  and  X is  the  wavelength  of  the  photons  in  A.  In  plotting  these 
expressions  in  Fig.  1,  they  were  converted  into  the  forms 


k(E) 


1.13  x 10  E 


-3.31 


(5) 


for  30  eV<E<282  eV  and 


k(E)  = 2.21  x 106E-3’68  (6) 

g 

for  282  eV<E<1550  eV,  using  the  definition  y = 10  (Ark/A).  The  expressions 

(3)  and  (4)  were  deduced  by  Lukirskii  et  al.  from  their  experimental 

2 

absorption  measurements.  The  results  of  Nordfors3  obtained  in  the  region 
from  1490  to  8050  eV  were  found3^  to  fit  well  to  the  expressions  (4)  or 
(6)  for  the  higher  energy  region.  Since,  as  will  be  shown,  in  the  spectral 
region  of  their  studies  the  refractive  indices  are  quite  close  to  unity, 
transmission  measurements  gave  the  absorption  coefficient  directly  without 
taking  into  account  reflection  at  the  film  boundaries.  Their  measurements 
were,  however,  made  at  only  eight  wavelengths  from  30  to  282  eV,  nine  wave- 
lengths from  282  to  1550  eV  and  f ive  wavelengths3-  from  1490  to  S050  eV, 
employing  the  characteristic  X-rays  from  various  anode  materials. 
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The  jump  in  absorption  it  282  eV  is  due  to  the  carbon  K-shell 
absorption  edge,  in  the  vicinity  of  which  some  fine  structure  should 
be  found.  The  study  of  fine  structure  in  core  shell  absorption  edges 
using  a continuum  source  of  synchrotron  radiation  is  of  current  interest  , 
but  such  studies  are  not  available  for  polystyrene.  Despite  the  coarseness 
of  Lukirskii  et  al.'s  and  Nordfors'  data,  their  results  seem  to  represent 
the  essential  features  and  correct  order  of  magnitude  of  the  absorption 
for  polystyrene  over  a very  wide  range  of  photon  energies.  Indeed  a 
fairly  close  fit  is  found  between  the  present  results  and  Lukirskii  et 
al.'s  results  in  the  low  energy  region.  The  difference  in  k between  the 
present  result  and  their  result  is  about  9%  at  30  eV,  and  that  between 
the  extrapolation  of  the  present  result  from  20  to  60  eV,  which  is  shown 
by  a dotted  line  in  Fig.  1,  and  their  result  is  about  8%  at  282  eV.  Two 
lines  representing  the  extrapolations  of  the  present  results  and  their 

I ^ 

results  Cwlv.O  eV  cross  at  around  100  eV,  indicating  that  the 
differences  are  not  due  to  any  systematic  errors  resulting  from,  for 
example,  error  in  determination  of  the  film  thickness  in  either  study, 
but  due  to  random  errors.  Their  original  data  at  eight  wavelengths 
below  282  eV,  in  fact,  exhibit  noticeable  scatter  from  a straight  line  fit 
even  on  a logarithmic  plot,  but  th  discrepancy  of  their  expression  (5) 
from  our  expression  (2.)  is  well  within  the  107  experimental  uncertainty 
c ia  i*  oc  by  Irte.t  or  thvir  i d ita. 

The  results  tor  k given  in  leg.  1 provide  basic  data  for  comprehensive 
infor:  • it  ion  about  the  electronic  response  of  polystyrene  to  electromagnetic 
radiation.  lie  fore  deriving  the  other  optical  quantities  from  these  results, 
the  validity  o!  th>-  !.  value  . obtained  hero  was  checked  by  evaluating  the 
su.:  rule* 
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N (E)  - • (m/4  n f.  -fi  ) f k ( I : ' ) d E ' >N 

L n o 


(/) 


where  n is  the  molecular  density  .m  a , i the  total  number  of  electrons 
in  a molecule.  Taking  a monomer  i i u i:  • of  polystyrene  as  a 

O o 

molecular  unit,  N i-  36.  In  carry  in  •.  euL  the  integration  (7)  , the  1 
values  above  S0;.n  eV  were  assumed  lo  be  represented  by  expression  (6)  up 
to  an  infinite  pi,  a on  energy.  Th i s can  be  expected  at  least  up  to  the 
r.  Lo ...  wi  i se  [ . : t -'mat  Lons  due  to  the  Compton  effect 

I 2 

and  ;>..  r-t  t i • ...pprec  iabln.  In  the  region  below  0.6  eV , where 

photoabsoi  ( : i ■ .Ocular  vibrations  occur  • , the  k values  were  assumed 

to  '.a  so  small  taat  tiie  contribution  of  E'k(K')  in  this  low  energy  region, 

to  N|(f;  could  be  neglected. 

ilculation  N.  ire  plotted  in  Fig.  2 as  a function 

of  E.  Since  in  this  study  it  was  difficult  to  determine  experimentally  the 

density  of  the  polystyrene  film,  p,  with  an  uncertainty  less  that  f-2/i,  the 

molecular  density  n^(=-L.p/M)  was  adjusted  so  that  the  value  of  N (tl)  became 

equal  to  N in  tfie  limit  of  K*"3.  (L  is  Avogadro's  number  and  M is  the 

3 

molecular  weight).  ThLs  gave  a film  density  of  1.040  g/cm  , which  is  in 

3 

excel:  .ut  agreement  with  the  measured  value  of  (1. 03±0.02)  g/cm  . Thus,  in 
this  check,  the  over-all  validity  of  the  k values  was  confirmed  with  an 
uncertainity  of  l/.. 

El  I t!wr  c ..i  lu  It  i'.cis  of  h (;.)  w-.  tin.  1c  on  the  k v.  ■ l • . du  » ' 1 t 

vale  . • and  core  exc  it  at  i >:i  • , epar.it  c ; y . ,\  saw-tooth  like'  shape  ci  the 
k »p  t rum  in  the  logar  ithmi.  p!  it  at  li  i ;h  etiorgv  a I 1 ows  us  to  separate 
t h<  I'  .arption  due  to  the  >o.  • It  of  i.  lie:"  the*  tot  i 1 absorption.  'ifie 
• " 1 ' ■ ' ] a : 1 > ; i v.  ; . • I «■  t 1 ! . . • i . i • : . 1 he  I e . t • 


9. 


K-shell  edge  to  higher  energies.  This  seams  to  be  a quite  plausible 

v c 

procedure.  Ttu  results  ot  Nj(E)  and  h'^(E),  the  contributions  to  Nj(E) 
Item  l h<*  v.ih'ui  i'  and  cat  e exr  itat  i our* , rouperl  ivelv*  are  given  in  Fig.  2. 
1 1 is  noted  that,  a redistribution  u!  electron  numbers  cm  responding 


to  1.51  electrons  is  found  between  the  valence  and  core  excitations, 
v c 

N^(E)  and  N^(E)  should  tend  to  40  and  16,  the  numbers  of  valence  and  core 

electrons,  respectively,  in  the  limitE  -*•  “>,  if  there  is  no  oscillator 

strength  coupling  between  them.  The  fact  that  the  oscillator  strengths 

for  inner  shells  are  considerably  smaller  than  their  numbers  of  electrons 

13 

has  been  known  since  the  early  days  of  X-ray  spectroscopy.  In  addition 
the  magnitude  of  the  oscillator  strength  redistribution  between  the  outer 
and  inner  shells  has  been  known,  in  general,  to  become  larger  as  the  atomic 
number  becomes  bigger.  It  is  further  to  be  noted  that  the  effective  numbers 
of  electrons  as  plotted  in  Fig.  2 as  functions  of  photon  energy  E do  not 
saturate  until  a few  thousand  eV.  This  is  so  even  for  the  valence  electrons 
showing  that  even  at  these  high  energies  the  contribution  of  valence  excitations 
to  the  oscillator  strength  is  not  negligible. 

The  analyses  presented  above  were  performed  on  the  present  results 
between  0.6  and  60  eV.  In  the  region  from  60  to  282  eV,  the  extrapolation 
of  the  present  results  between  20  and  60  eV  expressed  by  Eq.  (2)  was  used, 
instead  of  Lukirskii  et  al. 's  result  given  by  Eq.  (5)  for  the  reason 
mentioned  earlier.  Expression  (2)  was  also  used  to  represent  the  absorption 
by  v.-ticnce  electrons  above  23'  o'.  . I.  used  hr.k  Iraki  L »*•  al.  ’s  expression 
(5)  for  the  absorption  by  valence  electrons  above  60  el , no  substantial 
change  is  found  in  the  results  for  the  film  density  and  the  oscillator 
strength  redistribution.  In  all  further  analyses  presented  below,  therefore 
our  expression  (2)  was  used  for  the  absorption  by  v -vk,r\ce  electrons  above  60  eV. 


ID 


Some  uiu  ert  .1  i;U  u ■ . t lie 


i.  . .1  > 11,  iy  iri;e  from  the  ambiguity 


of  t ho  k data  i he 


a < • r p t i 1 cage 


at  around  282  eV . In  the  • I ■,  . edge  was  assumed 


at  282  eV  as  given  by  1 a. it.  . .in.  ■ sued  sharp  K-edge  has 


been  observed  for  polypropylene  > H •-),  an  r panic  polymer  similar 


to  polystyrene,  at  283  e\  im  soft  X ray  photons  front 


a synchrotron,  the  threshold  ener  ed  in  the  present  analyses  seems 


to  be  reasonable.  A ‘2  eV  deviation  the  threshold  t nergy  from  282  eV 


results  in  changes  1 nly  : 0.003  p/csn  and  +0. 13  elf.  trons  in  the  results 


for  the  : icillat  trength  redistribution,  respectively. 


Thus,  although  to  evaluate  the  uncertainty  involved  in  the  present  result; 


we  have  to  wait  for  a precise  observation  of  the  carbon  K-edge  structure  in  poly- 


styrene the  present  results  seem  to  be  sufficiently  c : rect  with  uncertainties 


pmtxxbh^ 

which  are  boss  than  those  suggested. 


2.  The  refractive  index  n. 


The  refractive  index  n(F.)  at  a ph  ' >n  energy  I'  can  be  determined  by  the 


Kramers-Kconig  relation 


»(r.)  - 1 


f E'k(E^) 


1 e prese 


. . , 


the  regions  out  tide  these  ierg|.  , the  values  v ' rnpolated  as 


■ 


described  earlier  in  the  sum  rule  calculation  can  be  used.  n(E),  therefore, 

can  be  determined  with  little  uncertainty  at  all  energies.  Just  as  k can 

be  separated  into  the  valence  and  core  contributions,  n(E)-l  is  also 

V c 

separable.  The  results  of  n (E)-l  and  n (E)-l,  the  contributions  to  n(F.)-l 
from  the  valence  and  core  excitations,  respectively,  are  presented  in  Figs. 

3 and  4,  respectively. 

V c 

Both  n (E)-l  and  n (E)-l  become  almost  constant  in  the  regions  far 

below  the  absorption  onsets.  It  is  to  be  noted  that  the  values  of  nC(E)-l 

below  about  100  eV  are  smaller  than  those  of  |nV(E)-l|  in  that  region  by 

3 

a factor  10  or  so,  showing  that  the  contribution  of  the  core  excitation 
to  the  refractive  index  is  negligibly  small  in  this  low  energy  region.  The 
singularity  appearing  in  n (E)-l  at  the  carbon  K-edge  is  due  to  the  present 

approximation  of  a step-like  onset  for  the  K-shell  absorption.  At  high 

V c 

energies,  n (E)-l  and  n (E)-l  tend  to  zero. 

As  will  be  shown  later  in  the  Appendix,  if  the  k(E)  values  are 
expressed  for  F.  >_  E^  by  a form 

k (E)  = ctF.-6  (9) 

as  has  indeed  been  found  in  the  present  results  for  E -20  eV,  n(E)-l  for 
K --  lOEj  can  be  written  as 

n (C)-l  = -Ej72E2  - Q(3)k(E)  (10) 

a 

with  an  accuracy  much  better  than  99/.,  where  F.  " is  a constant  defined  by 

P 

the  sum  rule  (7)  ns 
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E'k(E')dE' 


E 

P 


= (4 


2 ? 

n f I IN 
o l 


C 1 1) 


and  Q(8)  is  a qu  mt  Lty  which  depe  tl  ■ t !.  in  the 

expression  (9).  I t n in  tl  iix,  Figure  |£. 

shows  that  for  6 -2  as  in  ti  present  r.  _•  .Its,  n(E)--l  , ..es  asymptotically 
to 

n(E) -1  = -E  2/2V  (12) 

P 


at  high  energies.  This  form  of  n(  )-l  i t aat  of  a f i ee-electron  gas  in 

v . 

the  high  energy  limit.  In  the  . tl  e ..ions  k (E)= 


a _T  93  c 

7.81  x 10  E and  k (E) 

r 

n (E ) — 1 were  obtained  for 


- . : ' 10  '•  3-2'\  nV(F.)-l  and 

2820  t'V,  respectively,  as 


nV(E)-l  (i  V)-'  • -1.73  x 10'Vi  '’2\  (13) 

P 

'V"  • i • 78  x 103/E3’23, 

m 2 22  c 2 22. 

where  • > ■*'  I (4nn  e 1 /ra)  x 14.49.  The 

p o P o 

detailed  !■  *h  . 1 • : 1 t high  oner  . . a relative  to  the 

v , t c (•  7 

- r)  /2 

l' 

c wl  L i vt*  L j « .lit  1 > ) 1 • • . i.  ..it*  it.  t i i ,’S . 3 and  * * 

Ti  10.  t •.  I • • i - - U ! (a  (tj-l)/(r  (E)  -1 ) given  her 

r ' r 

nuitc  * s i»  i ; 1 1 ••  . < inver  ••  of  nV  (!'.)- 1 


and  900  eV,  respectively.  It  is  to  be  noted,  however,  that  the  deviation 
of  n (E)-l  values  from  the  f ree-electron  values  is  considerably  larger  than 
that  of  nV(E)-l  over  the  entire  high  energy  region.  The  relatively  close 
fit  of  the  nV(E)-l  values  to  the  f ree-electron  values  above  40  eV  in  Fig.  3 
has  to  be  considered  fortuitous.  As  is  shown  in  the  Appendix,  Figure  12, 
this  close  fit  is  due  to  the  relatively  small  value  of  Q(£5)  in  Eq.  (10)  for 
(5=3.23,  the  exponent  for  absorption  of  the  valence  electrons. 

A simple  and  interesting  sum  rule 


E 

S(E)  = { n( E ' ) -l}dE'->-0;  E -*•  °° 

A 


has  been  derived  recently  by  Altarelli  et  al.  for  the  refractive  index 
n(E)  o-f  isotropic  media  including  conductors.  This  sum  rule  asserts  that 
the  average  value  of  the  refractive  index  over  the  whole  spectral  range 
is  equal  to  unity.  In  the  derivation  of  this  sum  rule,  the  only  assumption 
made  was  that  in  the  high  energy  limit  the  medium  responds  like  a free- 
electron  gas,  i.e. , 


n(E)  + ik(E)  - 1-E  2/2E2. 

P 


As  we  have  already  seen,  this  is  indeed  the  case  for  the  present  results. 

V c 

S(E)  evaluated  from  n (E) — 1 and  n (E)-l  arc  given  in  Fig.  5 as  functions  of 

4 , <v , o 

photo  ■ energy  d up  to  10  eV.  Both  S ( .)  and  S (F.)  lonvi-rge  s’  ro  /.<.  rv 

at  high  energies.  The  analytical  expressions  (11)  and  (14)  allow  us  to 

v c 

extend  the  integrations  S (E)  and  S (E)  up  to  infinity.  It  was  found  that 

V z c 

both  S (E)  and  S (K)  converge  exactly  to  zero  in  the  limit  of  E 


i 
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For  a high  energy  E,  above  which  the  behavior  of  n(E)-l  can  bo 
represented  by  the  free-electron  expression  (12),  S(E)  may  be  approximated 
by 


S (E) 
F 


(17) 


This  form  describes  the  asymptotic  behavior  of  S(E)  in  the  high  energy  limit. 
Sp(E)  from  the  valence  and  core  contributions  are  presented  in  Fig.  5 and 
compared  with  the  results  for  S(E).  S(E)  is  described  well  by  S (K)  over  a 
wide  range  in  the  high  energy  region.  Relative  deviations  between  S(F)  and 
S^,(E)  are,  however,  about  the  same  as  those  found  in  the  results  for  (n(E)-l)/ 
(nf.(E)-l^at  corresponding  energies  as  shown  in  the  inserts  of  Figs.  3 and  4. 

Before  closing  this  subsection,  we  present  in  Fig.  6 the  results  of 
n(E)-l,  in  the  region  where  n(E)-l  is  negative,  on  logarithmic  plots, 
together  with  the  free-electron  values.  These  plots  are  of  help  in  evaluating 
the  complex  dielectric  function  and  the  energy  loss  function,  which  will  be 
presented  in  the  following  section. 


3.  The  complex  dielectric  function  e and  the  energy  loss  function  Ir.(-l/e). 

Toe  imaginary  part  of  the  complex  dielectric,  function  *(-••- : +i:,>  r. r 
energy  loss  function  Im(-l/e)  for  fast  charged  particles  were  calculated 
from  the  values  of  n and  k,  and  are  presented  in  Fig.  7 on  logarithmic  plots, 
together  with  the  values  of  2k.  Among  these  three  spectra,  substantial 


differences  are  found  only  in  the  low  energy  region.  Deviations  of  the 
Imfl/e)  spectrum  from  » found  in  tiie  low  energy  region,  describe  the 
so-called  collective  effect,  J i.e.  charged  particle  density  fluctuations 


brought  about  by  fast  charged  particles  incident  on  condensed  natter.  In 


regularities  in  the  displacement  of  the  Low  energy  structures  should  be 
noted.  All  three  peaks  found  in  the  low  energy  region  shift  to  higher 
energies  in  order  of  £ 9 , 2k  and  Im(-l/r),  and  their  peak  heights  decrease 
also  in  this  order. 

For  and  Imfl/t),  there  exist  the  well-known  f-sum  ruLes,^" 


n2(£) 


(m/ATin^e2^2)  j E'e2  (E')dE'->-N;  E-*», 
o 


(i«> 


N3(E) 


(m/4TTn 
tt  o 


fE 

E'lm  {-1/ e (E')  }dE'->N  ; 
o 


E-*», 


(l^) 


The  results  of  N2(E)  and  N^(E)  are  shown  in  Fig.  8 as  functions  of 
photon  energy  E,  together  with  that  of  N (E)  which  has  already  been 
presented  in  Fig.  2.  It  was  found  that  all  these  numbers  of  effective 
electrons  converge  to  the  same  value  in  the  high  energy  limit.  As  has 
been  seen,  k and  n-1  are  separable  into  the  valence  and  core  contributions. 
Generally,  e2  and  Im(-l/e)  are  not  separable,  as  by  definition  they 
include  the  cross  terms  between  the  valence  and  core  contributions.  In 
the  case  of  |n  -1|  and  | n —1 | << 1 , and  k and  k <<1,  however,  it  can  be 
show  that  f9  and  Im(-l/r)  are  separable.  This  can  be  seen  from  the 
results  shown  in  Fig.  7,  where  f and  Im(-l/t)  are  in  agreement  with  2k 
in  the  high  energy  region;  thus  they  can  be  separated  into  the  valence 
and  cot  contr  i but  Inns  in  the  same  manner  as  th  1 spoei  run.  i.  ( (!.  ■ , 

N2(F)  and  N((E)  fur  the  valence  and  core  contributions  are  plotted  in 
Fig.  8.  because  of  the  close  agreement  among  2k,  t 2 and  lm(  • 1 /•  ) , no 
substantial  differences,  were  found  in  the  N(E)  values  in  t Ite-  high  energy 
region. 

Lot  * 
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Detailed  spectra  for  i j,  » } , ami  1 in (— 1 / < ) in  the  low  energy  region 
are  presented  in  Figs.  l)  and  10  on  linear  plots.  The  general  spectral 
features  given  here  agree  well  with  results  from  earlier  studies'^  ’ ^ ^ 
obtained  below  ?0  eV.  Considerable  numerical  disagreements  are,  however, 
found  between  the  present  results  and  some  of  the  earlier  ones.  Apparently, 
the  discrepancies  can  be  attributed  to  poor  experimental  accuracy  and  large 

uncertainties  in  data  analysis  in  the  earlier  studies.  Four  absorption 

bands  can  be  distinguished  below  11  eV  in  the  spectrum  shown  in  Fig.  9; 

at  around  4.8,  5.8,  6.3,  and  9.5  eV.  The  interpretation  of  these  structures 

and  the  corresponding  structures  in  Im(-l/e)  have  already  been  described 
9 16-20 

by  many  authors.  ’ The  spectral  features  in  the  extreme  ultraviolet 

above  10  eV  given  in  Fig.  10  are  quite  simple  and  are  marked  only  by  the 

broad  peaks  at  12.5  eV  in  and  at  21.7  eV  in  Im(-l/e).  These  features 

6 7 21  22 

are  typical  of  many  organic  solids.  ’ ’ ’ The  spectrum  for  Im(-l/e) 

shown  here  agrees  well  with  that  from  characteristic  electron  energy  loss 
measurements^ ’ ^ on  thin  polystyrene  films. 

4.  The  average  photoexcitation  energy  and  the  average  energy  loss. 

The  present  results,  extending  substantially  over  the  region  from 
zero  to  infinite  photon  energy,  allow  us  to  evaluate  the  average  photo- 
excitation energy  and  the  average  energy  loss  suffered  by  fast  charged 
particles.  These  quantitites  are  of  key  importance  in  radiation  physics 
which  deals  with  the  interaction  of  high  energy  ionizing  radiations  with 
matter . 

Shown  in  Fig.  11  are  the  results  of  the  photoexcitation  energy 
£2^^)  and  the  energy  loss  F.^(K^)  averaged  over  the  spectral  region  from 
Kq  to  infinity,  i.e.. 
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W = fi2(E)dE/  c,(E)dE, 


E3(Eo)  - E Im{ -1/t (E)  }d£/  Ira{-l/t.(E)  }dE.  (2() 


Such  plots  ot  E (E  ) and  E (E  ) as  functions  of  E are  of  help  in  evaluating 
2 o 3 o o 6 

the  contribution  of  excitations  in  a particular  energy  range  to  the  average 
excitation  energy  and  energy  loss.  Since,  in  the  present  results,  the  onset 
of  valence  absorption  was  found  at  4.5  eV,  the  values  of  E^(E{))  and  E^CE^) 
at  E^  = 4.5  eV  give  the  average  excitation  energy  and  energy  loss, 
respectively,  over  the  entire  oscillator  strength  distribution.  E2  and  F.^ 

at  E^  = 4.5  eV  were  found  to  be  25.1  and  36.8  eV,  respectively.  This 
difference  can  be  attributed  to  an  energy  loss  associated  with  the  creation 
of  collective  oscillations  due  to  the  interaction  of  fast  charged  particles 
with  matter.  Photons  generally  do  not  excite  collective  oscillations  in 
experimental  systems  designed  for  the  measurement  of  k. 

From  the  definition  of  the  energy  loss  function,^’  the  average  energy 
loss  F-3(E^)  obtained  here  should  be  considered  to  be  that  due  to  a single 
energy  loss  event  associated  with  zero  transverse  momentum-transfer  from 
the  fast  charged  particle  to  the  medium.  Therefore,  it  cannot  be  compared 


direct1 y with  similar  quantities  such  as  the  mean  excitat 


ton  potent  tal  , 


24  '>  ^ 

(O  , appear  ing  in  the  wll-known  ilethe-Uloch  tonnula  ’ tor  the  .. capping 
power  or  the  so-called  W-value,  the  average  energy  for  an  ion-pa iv  creation 


It  may  be  worthwhile,  however,  to  note  that  f = 63.5  <>V  was  obtained  hv 


Sternheimer  for  polystyrene,  and  W = 37  oV  was  reported  by  Mathieu 
27 

et . al.  for  liquid  cyclohexane,  which  may  be  considered  to  be  a 


model  system  for  solid  polystyrene. 

The  variations  of  E2(Eq)  and  E3(Eq)  shown  in  Fig.  11  as  functions 

of  E are  no  more  than  a few  eV  for  E < 10  eV,  showing  that  the 
o o - 

contributions  of  such  low  energy  excitations  to  the  average  excitation 

energy  and  energy  loss  are  relatively  small.  This  result  will  apply  to 

other  condensed  hydrocarbons,  as  substantial  differences  of  the  oscillator 

21 

strength  distribution  for  these  materials  have  been  found  only  in  the 
region  below  about  10  eV.  Also  shown  in  Fig.  11  is  E^(Eq)  defined  by 


E (F.q)  = Ek(E)dE  / k(E)dE. 


As  is  well-known,  Ek(E)  for  condensed  media  cannot  be  considered  to  give 

the  oscillator  strength  distribution  in  a strict  sense.  It  has  been 

9 28  29 

found  for  many  hydrocarbons,  * ’ however,  that  except  for  molecular 

Rydberg  fine  structure  the  k spectrum  in  the  gas  phase  agrees  well  with 

that  in  the  condensed  phase,  apart  from  a constant  factor  corresponding 

to  the  ratio  of  molecular  densities  between  these  two  phases.  This  is 

particularly  so  in  the  region  above  10  eV,  where  most  of  the  excitation 

energy  and  energy  Loss  occurs.  !'  (E  ) shown  her-,  therefore,  can  he 

1 o 

regarded  as  representing  approximate  values  of  E„(E  ) and  F,„(F.  ) in  the 

l o do 

low  density  limit  i.e.  for  gas  phase  hydrocarbons,  since  f and  lm(-l/c) 


j i > 
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in  the  low  density  limit  agree  with  2k  over  the  whole  spectral  range^p 

^ E^V(Eo>,  F^CE^)  and  E^  (F.^)  presented  in  Fig.  11 
are  the  respective  average  excitation  energy  and  energy  loss  derived 
from  the  valence  excitation  only.  Here,  one  may  see  how  much  the  core 
excitations  contribute  to  the  average  excitation  energy  and  energy  loss. 

Finally  the  values  of  the  extinction  coefficient  k obtained  and 
used  in  this  study  are  tabulated  in  Table  I for  further  use. 


5.  Concluding  remarks. 

It  was  found  in  the  present  study  that  the  convergence  of  all  sum 
rules  associated  with  the  optical  constants  of  polystyrene  are  quite 
slow,  and  the  almost  perfect  saturations  do  not  occur  until  the  keV  region. 
This  is  for  a material  having  one  of  the  simplest  of  all  electronic 
configurations.  Experimentally,  the  f-sum  rules  have  provided  useful 
criteria  for  consistency  checks  of  the  optical  constants.  It  has  been 
shown,  however,  that  the  f-sum  rules  are  applicable  only  when  the 
experimental  data  cover  all  electronic  excitations  including  the  innermost 
core  excitations.  In  applying  the  f-sura  rules  for  a particular  shell, 
difficulties  arise  from,  in  addition  to  the  slow  convergences,  the  fact 
that  a non-negligible  amount  of  oscillator  strength  coupling  always  occurs 
between  the  inner  and  outer  shell  electrons.  If  experimental  data  cover 
only  a part  of  the  total  oscillator  strength  distribution  as  is  so  ■>  r. 
most  cases,  the  oscillator  strength  for  a particular  shell,  to  which  the 
f-sum  rules  should  converge,  can  not  be  known  a priori  without  knowing 
those  for  the  remaining  shells. 
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It  h.is  been  proposed  recently  by  AU.irolli  ot  al.w  th.it  the  S(F.)  sum 
rule  for  the  refractive  index  given  by  Fq.  (15)  may  be  used  as  a 
saturation  criterion  and/or  as  a consistency  check  for  the  optical 
constants  obtained  experimentally.  The  present  results  for  S(E)  shown 

in  Fig.  5,  indeed,  converge  to  the  f ree-electron  values  S (E)  given 

r 

by  F.q.  (17)  much  faster  than  the  f-sum  rules  converge.  Differences 
between  S(E)  and  S (E)  , however,  are  not  negligibly  small  up  to  quite 

r 

2 

high  energies.  Further,  E^  for  a particular  shell  in  the  expression 
for  Sp(E) , which  is  proportional  to  the  oscillator  strength,  can  not 
be  know  until  the  convergence  of  the  f-sum  rule  is  attained.  Practically, 
the  new  sum  rule  for  S(E),  if  applied  to  a particular  shell,  is  subject  to 
the  same  uncertainty  as  the  f-sum  rules. 


APPENDIX 

Let  us  assume  that  the  extinction  coefficient  k(E)  at  photon  energy 
E is  expressed  by 

k(E)  = aF."6  (Al) 

for  E>E  as  is  generally  found  for  high  photon  energies,  where  a and 
P are  positive  constants.  According  to  a Kramers-Kronig  relation,  the 
ret:  t i v.  • Index  n(h)  at  photon  rgy  K is  givo.u  bv 


n(F.)  -1 


2 


T. 


E ' k ( E * ) 
. F - E2 


dF. 


(A2) 
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n. 

Dividing  the  integration  into  three  regions,  and,  then,  allowing  less;  than 
1%  error  to  the  first  and  third  regions,  n(E)-l  can  be  written  as 


n (E) -1  = — - 

IT 


f 0.  IE 
JQ 


E'k(t: ') 
E'z-E7 


dE'  + 


rlOE 
0 . IE 


E'k(E') 

e 


dE 


+ 


E'k(E') 


10E 


-E 


dE' 


; 


0.  IE 

E'k(E')dE'  + 

0 


rlOE 
0.  IE 


E'k(E') 

E7y~F/ 


dE' 


k(E') 

E' 


dE'}. 


(A3) 


Now,  if  we  confine  ourselves  to  the  case  of  E>10E  . The  first  integral  can 
be  written  as 


~ (l/E2)  {- 


r “ 


E'k(E')d£'  + 


E'k(E')dE'). 


(A4) 


0.  IE 


The  first  integral  in  (A4)  should  be  a constant  from  the  sum  rule  , and  the 
second  one  is  calculable  analytically,  as  the  expression  (Al)  holds  for 
this  region.  Thus,  the  integrals  in  (A4)  become 


-E2/2E 

P 


2 


10 


8-2 


6-2 


k(E)  , 


(A5) 


where  E is  defined  by 
P 


E'k(F.')dE' 

o 


(A6) 


: d integral 


in  Eq . (A3)  is  also  calculable  analytically  and  i.-.  giver. 


2 

t; 


. 

(A7) 


01 


Finally,  the  second  integral  in  Kq.  (A3)  is  also  calculable  as  the 
expression  (Al)  held,  foi  t 1 1 i : . legien.  lloWt'VOl,  the  luu.tienil 
dependence  ol  this  integral  on  l.  can  be  obtained  easily  without 
carrying  out  the  integration.  Using  the  similarity  o.  the  expression 
(Al)  in  logarithmic  scale,  the  integral  region  from  0.1E  to  10E  can 
be  shifted  to  an  arbitrary  region,  for  example,  from  O.IE^  to  10E2. 
This  shift  introduces  a factor  (E/E  ) i.e.. 


7 r10E 


- J0.1E 


E'k(E') 
E ^ -Et 


dE' 


(e/e2)  b 


r10E2 

■'o.ie2 


E'k(E') 


dE' 


(A8) 


Since  this  shift  of  integral  region  is  valid  for  any  E-2(>E),  the  integral 


(i/e2)  b 


dE 


(A9) 


should  be  a constant.  Thus,  the  second  integral  in  Eq . (A3)  can  be 
written  as 

-~r  P(S)k(E)  (A10) 


where  ?(S)  defined  by 

o f10Ey  ,.^-E+l 

P(0)  = O/E  ) I E (A  LI ) 

JQ.1E2  ~l‘2 


i a unit  i t •,  de  ••aa  i i,,  oa  i v on  ...  Su.nrtuig  up  (,.\5)  , (A7),  and  (A  10  J , we 

have  for  E>10-;^ 


i 

t 

i 


1 


! 


H 


n ( E ) - 1 = -(t:  ?/2E?‘)+  — 
P '• 


, 

6-2 


+ P(8) 


10 

6 


-6 


ik(F.) 


(AL2) 


Values  of  Q(b)  obtained  from  a numerical  calculation  of  I'(fS)  are  plotted 
in  F ig . \z  as  a function  or  6.  Q(S)  becomes  zero  at  S-3.06  and  goes  to 
negative  infinity  in  the  limit  of  g>2.  Such  an  analytical  expression  for 
n(E)-l  for  the  high  energy  region  is  quite  useful  in  analyses  of  optical 
data  as  has  been  show  in  this  study,  and  can  he  applicable  for  any  kind 
of  material,  as  the  expression  (Al)  involves  a general  behavior  of  k due 
to  high  energy  absorptions. 
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FIGURE  CAPTIONS 


Fig.  1. 


Fig.  2. 


Fig.  3. 


Fig.  4. 


Fig.  5. 


The  extinction  coefficient  k for  polystyrene  as  a function  of 
photon  energy  E.  Dotted  line  is  extrapolation  of  the  present 
results  between  20  and  60  eV. 

Effective  number  of  electrons  N^(E)  per  monomeric  unit  of 

polystyrene  obtained  from  the  extinction  coefficient  k as  a 

v c 

function  of  photon  energy  E.  N^(E)  and  N^CE)  represent  the 
contributions  from  the  valence  and  core  electron  excitations, 
respectively . 

The  refractive  index  for  polystyrene  due  to  the  valence  electron 
excitations  as  a function  of  photon  energy  E.  Dashed  curve 
represents  the  f ree-elec tron  values. 

The  refractive  index  for  polystyrene  due  to  the  carbon  core 
electron  excitations  as  a function  of  photon  energy  E.  Dashed 
curved  represents  the  free-electron  values. 

The  sum  rule  S(E)  for  the  refractive  index  n of  polystyrene  as 

a function  of  photon  energy  E.  SV(E)  and  SC(F.)  represent  the 

contributions  from  the  valence  and  core  electron  excitations, 

V c 

respectively.  S (E)  and  S (E)  represent  the  corresponding 
r r 

free-electron  values. 

The  refractive  index  n for  polystyrene  in  the  region  where  n-1 

is  negative  as  a function  of  photon  energy  E.  nV  and  nC  represent 

the  contributions  from  the  valence  and  core  electron  excitations, 

V c 

respectively,  n^,,  n^,  and  represent  the  corresponding  free- 
electron  values. 
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Fig.  7. 


Fig.  8. 


Fig.  9. 


Fig.  10. 


Fig.  11. 


Fig.  12. 


Twice  the  extinction  coefficient  k,  the  imaginary  part  of  the 
complex  dielectric  function  t.  (-£  y + it  ) and  the  energy  loss 
function  Tm(-]/r)  for  polystyrene  as  functions  of  photon  energy  F.. 
Effective  number  of  electrons  N ^ (E) , N2(E),  and  N3(E)  per  monomeric 
unit  of  polysytrene  obtained  from  the  extinction  coefficient  k,  the 
imaginary  part  of  the  complex  dielectric  function  e^Ej+ie^)  and 
the  energy  loss  function  Im(-l/e),  respectively,  as  functions  of 
photon  energy  E.  NV(E)  and  NC(E)  represent  the  contributions 
from  the  valence  and  core  electron  excitations,  respectively. 

The  real  and  imaginary  parts  of  the  complex  dielectric  function 
£(=£^+^2)  and  the  energy  loss  function  Im(-l/e)  for  polysytrene 
as  functions  of  photon  energy  E. 

The  real  and  imaginary  parts  of  the  complex  dielectric  function 
£(=Ej+it2)  and  the  energy  loss  function  Im(-l/e)  for  polysytrene 
as  functions  of  photon  energy  E. 

The  average  photoexcitation  energy  E2  and  the  average  energy  loss 
E^  for  polystyrene  over  the  spectral  range  from  EQ  to  infinity 
as  functions  of  E^.  may  represent  E2  and  E^  in  the  low  density 

— v — v — V 

limit.  E^,  E2>  and  E^  represent  the  corresponding  values  due  to 
the  valence  electron  excitations  only. 

Dependence  of  Q(S)  on  the  exponent  8. 
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Table  I.  The  extinction  coefficient  k for  polyst  y rt*ne 
energy  E . 


E (eV) 

k 

E(eV) 

k 

0.6 

11.0 

0.586 

: 

0.00000 

11.5 

0.636 

4.4 

12.0 

0.681 

A. 5 

0.0011 

12.5 

0.713 

A. 6 

0.0089 

13.0 

0.730 

A. 7 

0.0103 

13.5 

0.743 

A. 8 

0.0103 

14.0 

0.748 

A. 9 

0.0092 

14.5 

0.747 

5.0 

0.0072 

15.0 

0.741 

5.1 

0.005A 

16.0 

0.702 

5.2 

0.00A1 

17.0 

0.646 

5.3 

0.0043 

18.0 

0.590 

5. A 

0.0292 

19.0 

0.531 

5.6 

0. 183 

20.0 

0.473 

5.8 

0.330 

22.0 

0.361 

6.0 

0.A31 

24.0 

0.276 

6.2 

0.825 

26.0 

0.211 

6. A 

1.023 

28.0 

0.164 

6.6 

0.88A 

30.0 

0.129 

6.8 

0.694 

35.0 

0.0806 

7.0 

0. 524 

40.0 

0.0542 

7.3 

0.354 

45.0 

0.0375 

7.6 

0.239 

50.0 

0.0247 

7.9 

0.206 

55.0 

0.0192 

8.2 

0.220 

60.0 

0.0136 

8.6 

0.264 

: 

a 

9.0 

0.344 

282 

9.5 

0.425 

: 

b 

10.0 

0.474 

8050 

10.5 

0.523 

k=7.81xl03K~ 3 

2 3 

(Extrapolation  of 

k between  20 

and  60 

k=2.21xl06E~3 

(Refs.  2 and  3; . 
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APPENDIX  C, 

Reprinted  from:  Health  Physics  Division  Annual  Progress  Report 

for  Period  Ending  June  10,  197b,  ORNL-5171,  pp.  160-162. 

ELECTRON  ATTENUATION  LENGTHS  IN  THE 
VACUUM  ULTRAVIOLET  AND  SOFT 
X-RAY  REGIONS1 2 

Considerable  time  has  been  spent  this  year  reviewing  deposited  overlayers  is  subject  to  a practical  lower  limit 

attempts,  both  by  us1 1 and  by  others  as  reported  in  the  of  about  200  eV  for  the  electron  energies  available.  It 

literature,  to  obtain  electron  attenuation  lengths,  L,  appeared  that  hotoemission  was  a promising  approach 

from  observations  of  photoelectron  yield  in  the  ultra-  to  obtaining  /.  i , this  energy  range, 

violet  and  soft  x-ray  regions.  There  is  great  theoretical  Over  the  past  several  years,  we  have  attempted  to 

and  practical  interest  in  the  related  quantities,  the  mean  obtain  electron  attenuation  lengths  from  observations 

free  path  of  electrons  in  the  bulk  of  a solid  and  the  of  the  total  photoelectric  yield.  Y(0 ),  as  a function  of 

escape  depth,  range,  or  attenuation  length  for  electrons  the  photon  angle  of  incidence.  0.  Last  year13  we 

emitted  across  the  surface.  Observations  on  electrons  reported  that  a theoretical  analysis  showed  large  uncer- 

with  energies  from  a few  electron  volts  above  the  Fermi  tainties  in  the  L values  calculated  from  our  observation 

energy  to  about  200  eV  are  of  particular  interest,  as  of  G(0)  vs  0,  where  (7(0 ) = Y(6)!Y{ 0).  These  observa- 

lheones  predict  a variation  with  energy  of  the  electron  tions  were  for  thick  vapor-deposited  carbon  films  over 

mean  free  path  which  has  a minimum  in  the  region  of  the  range  of  photon  energies  from  20  to  64  eV.  where  L 

50  to  100  eV  for  mos*  materials  It  is  just  this  energy  is  known  to  be  small  In  fact,  for  any  material,  if  l.  is 

region  which  has  been  least  studied  and  where,  conse-  about  10  A or  less  as  is  generally  thought  over  the 

quently,  there  are  relatively  few  published  values  for  range  of  electron  energies  from  about  10  eV  to  about 

electron  mean  free  paths  in  solids.  In  the  ultraviolet  and  1000  eV,  then  to  obtain  L to  ± 1 A,  G(0),  the 

soft  x-ray  regions,  materials  are  either  strongly  or  extinction  coefficient  k,  and  the  polarization  P of  the 

moderately  absorbing  over  at  least  part  of  the  energy  incident  photons  must  each  be  known  to  about  1%.  For 

range.  At  the  same  time,  electron  attenuation  lengths,  our  observations,  both  G(0)  and  Fare  reproducible  and 

L,  are  short.  S 10  A,  over  much  of  the  region.  Measure-  can  be  obtained  to  ±1%.  However,  for  materials  wi 

ment  of  these  short  attenuation  lengths  requires  special  are  moderately  to  strongly  absorbing,  a limit  of  error  on 
techniques.  Two  methods  frequently  used  to  measure  k of  ±1%  is  unattainable  at  this  time.  In  attempting  to 

electron  attenuation  lengths  are  obviously  not  appli-  calculate  L from  observations  of  G(0)  vs 0,  inaccuracies 

cable  in  this  case  Electron  transmission  measurements  may  also  lesult  if  the  sample  is  not  homogeneous  with  a 

require  homogeneous  self-supporting  thin  films.  For  clean,  smooth,  planar  surface.  In  addition,  there  are  still 

very  short  attenuation  lengths,  it  would  seem  to  be  questions  concerning  the  theory  by  which  L is  related 

impossible  to  prepare  such  films  which  are  sufficiently  to  the  measured  G(0)  vs  0.  Different  theories  yield 

thiii.  The  use  of  Auger  electrons  generated  in  a significantly  different  L values  for  the  same  input 

substrate  to  determine  electron  attenuation  lengths  in  data.14  Finally,  since  the  definition  of  L assumes 

homogeneous  media  and  L may  be  considered  a 
• macroscopic  quantity  even  if  the  experimental  and 

12  Research  sponsored  in  pan  hvthe  Defense  Nuclear  Agency.  theoretical  problems  are  dealt  with  satisfactorily,  there 

under  Subtask  TA040.  but  does  not  reflect  endorsement  by  the  remains  the  philosophical  question  of  the  significance 

13  fc  i Arakawa  ei  al . Heahh  Phyn  D,v . Annu  ProK  ReP.  of  i,  values  that  are  calculated  to  be  comparable  with  or 

June  io,  1975.  OR\t.-5046,  pp  I25~2X.  and  previous  reports  m less  than  the  surface  roughness  of  the  sample. 

'bis  senes  It  was  thought  that  an  alternative  method  of  ob- 

u x i Khewer.  preprint  taining  L from  photoemission  measurements  would  be 

analysis  of  Y(t)  vs  the  film  thickness  r.  The  range  of 
film  thickness  required  depends  on  the  optical  proper- 
ties of  the  film  and  the  value  of  L of  the  emitted 
electrons.  Fot  a strongly  absorbing  film  and  small  L,  the 
thicknesses  available  should  range  from  t < L to  1 = <*>. 
For  weakly  absorbing  flints,  one  does  not  need  to  have 
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films  with  / v.  / . but  small  enough  thicknesses  should 
be  used  that  the  first  and  strongest,  structure  m Mr)  vs 
! is  accurately  specified.  Furthermore,  these  films  need 
to  be  homogeneous  and  continuous  and  to  have  smooth 
interlaces  with  the  substrate  and  with  vacuum  in  order 
for  MO  vs  t to  yield  meaningful  / values.  A literature 
survey  showed  that  usuig  currently  available  techniques 
for  making  thin  films,  all  films  are  discontinuous  for  i 
less  than  about  50  to  300  A.  the  thickness  at  which 
continuity  occurs  depending  on  the  method  of  prepara- 
tion. In  many  cases,  values  of  I.  calculated  from 
observations  on  ultrathin  films  would  be  more  de- 
pendent on  the  method  of  film  preparation,  and  hence 
film  structure,  than  on  the  basic  electronic  properties  of 
the  material  ol  the  film.  Thus  this  method  is  generally 
unsuitable  for  the  determination  of  I.  values  in  the  uv 
and  soli  x ray  region  In  view  ol  this  conclusion,  some 
/ values  obtained  bv  ibis  method  and  published  m the 
literature  turn  out  to  be  meaningless,  except  insolai  .is 
they  set  upper  bounds  on  the  aetual  values  of  / 

Examination  of  variations  ol  these  two  methods,  such 
as  (7(0)  vs  0 for  thin  films  of  known  thickness,  or  the 
ratio  of  photoelectric  yields  for  a thin  film  illuminated 
trorn  Ihe  Iront  and  Irom  die  back  through  the 
substrate,  as  a function  of  film  thickness,  reveals  diat  all 
methods  tor  /.  involving  photoelectric  yield  measure- 
ments appear  subject  to  the  same  limitations  when  L is 
small  and  the  material  is  moderately  to  strongly 
absorbing 

In  conclusion,  when  the  merirum  is  strongly  ab- 
sorbing. techniques  involving  thick  dims  must  be  used 
to  determine  I values.  However,  to  obtain  A to  H A 
when  /.  is  less  than  about  10  A requires  that  k be 
known  to  within  ±1'T.  For  strongly  absorbing  materials 
in  the  energy  range  of  interest,  k cannot  presently  be 
measured  to  this  accuracy.  Thus  L cannot  be  obtained 
with  any  meaningful  degree  of  accuracy  When  the 
medium  is  weakly  absoibtng  and  /.  is  short,  both  thick  - 
and  thin-fiim  techniques  can  he  used,  but  data  should 
be  interpreted  with  great  care  Films  should  not  be 
employed  that  are  so  thin  riut  they  are  not  continuous. 
I tfective  /.  values  calculated  from  observations  on  such 
films  are  dependent  on  film  structure  and  hence  are  not 
characteristic  of  the  material  When  effective  /.  values 
cati  be  obtained  by  a method  employing  thin  dims  that 
are  (luck  enough  to  be  continuous  and  homogeneous, 
they  should  tie  corrected  tot  si/e  ettects  so  that  the 
electron  attenuation  lengths  obtained  correspond  to 
those  lor  electrons  escaping  across  a planar  interlace 
between  die  semi-infinite  medium  and  vacuum.  Using 
currently  developed  experimental  techniques,  for 


weakly  absorbing  mateiial  I.  values  can  he  obtained, 
but  it  is  extremely  difficult  to  get  values  accurate  to 
even  t Iri  k by  any  photoeniission  technique.  Finally, 
especially  for  short  attenuation  lengths,  the  theories  of 
photoeniission  employed  in  their  derivation  need  to  be 
refined. 

PHOTOEMISSION  FROM  THICK  SAMPLES'  ‘ 

Although,  as  described  in  the  preceding  section, 
observations  of  (AO ) vs  0 from  a (luck  sample  cannot  be 
analyzed  for  electron  attenuation  lengths  when  the 
sample  is  strongly  absorbing,  the  relative  yields  for  a 
given  photon  energy  are  found  to  be  characteristic  of 
(lie  material.  Generally,  (7(0)  vs  0 is  reproducible  to 
tl%  and  is  relatively  insensitive  to  surface  conditions. 
These  data,  currently  unavailable  in  the  open  literature, 
are  of  practical  interest  in  the  development  of  improved 
plioloenntters  m the  vacuum  ultraviolet  and  soft  x-ray 
regions.  Over  the  past  several  years  we  have  obtained 
observations  of  (1(0)  vs  0 tor  insulators  such  as 
polystyrene  and  amorphous  carbon,  semiconductors 
such  as  silicon,  and  meials  such  as  gold  and  aluminum 
at  selected  energies  in  ihe  range  of  photon  energies 
from  about  20  to  about  SO  eV. 

Representative  data  for  carbon  have  been  presented 
previously.15  Representative  data  for  insulators,  semi- 
conductors, and  inetals  are  given  in  Figs.  29.3-29.5.  As 
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t ig.  29.3.  Relative  photovield  vs  angle  of  incidence  for  three 
different  energies  of  photons  on  polystyrene. 
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t ig.  29.4.  Relative  pholoyietd  vs  angle  oi  incidence  fot  three 
different  energies  of  photons  on  silicon. 


Ilian  a lacloi  it!  10,  ihe  magnitude  o!  O(0nl)  is  seen  lo 
lie  dommaled  by  the  value  ot  A 
In  view  of  the  conclusion  teaclied  ill  the  |ncvnius 
section,  namely,  dial  out  0(0)  vs  f / data  in  (lie  icguul 
liom  .’()  lo  KC)  eV  cannot  be  analyzed  at  the  piesenl 
time  lo  yield  meaningful  /.  values  all  of  our  accumu- 
lated 0(0)  vs  0 data  aie  currently  being  piepaied  for 
publication  without  attempting  to  obtain  the  associated 
/.  values. 
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Fig.  29.5.  Relative  pliotoyicld  vs  angle  of  incidence  for  39-cV 
photons  on  prld.  At  39  eV,  the  optical  functions  of  gold  ate  n = 
0 77  and  k 0 41.  The  polarization  was /'  - 0.87 

expected,13  the  characteristics  ot  0(0)  vs  0 are  de- 
pendent in  each  case  at  any  given  energy  on  the  optical 
functions,  n and  A,  of  the  material  at  that  energy  and 
the  polarization  of  the  incident  light  The  refractive 
index,  rt,  determines  the  position  of  the  maximum  in 
0(0)  by  the  relationship  n sin  0m,  while  the 
extinction  coefficient,  k,  and  electron  attenuation 
length,  A.  determine  the  magnitude  of  0(0, „ ).  Since  for 
all  the  data  presented  it  is  believed  that  /.  is  small  and 
for  any  practical  interpretation  cannot  vary  by  more 
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lAum  AIDIN'.  in  MI  AN  I Kit  i'AIhS  AND  STOPPING  POWERS  Of  ion  l NLRUY  1 1 1 1: 1 RONS 
ID  keV  > IN  SOLIDS  usm;  A STATIST  ICAL  MODI  L • 

' Ashley,  l lung  , R.  H Ritchie,  and  V.  f Anderson' 
th  i’hyst.s  > . '>iov  Oak  Ridge  National  laboratory 
Oak  Ridge,  Tennessee  37830 


Abstract 


II  Electron  Gas  tatistical  Model 


A statistical  inode  ' is  described  and  enip  1 ov » j 
calcu’ate  inverse  mean  tree  paths  and  stopping  powers 
fpr  electrons  of  energies  from  a few  eV  to  10  ce.  a; 
tne  pen*i'  Ipve)  in  \)  t Si,  Ni  , Cu,  Ad),  and  An.  b* 
tab  es  : f near  free  paths  and  stopping  powers  V» 
v i ids  d-T  presented  In  some  cases  ..apt-’  al  c i 
■ * nver  t-  '-.i-afi  f»*t»  path*,  and  stoppinn  powe  r*  are  a 
in  ' ude  ! ’he  calcolat iff.s  oased  on  this  mooei  arf* 
discussel  and  compared  *.-*‘th  previous  work 

troduction 


> e n era  ^ Theory 

A (.hanged  particle  passing  through  a solid  inter- 
a i with  a large  number  of  electrons  simultaneously 
and  it  is  thus  appropriate  to  speak  of  a mean  free  path 

• '*>p  charged  particle  against  enerqy  loss  to  the 

r 1 Assuming  the  effe<  t of  the  charqed  particle  on 
••  r^diun.  may  t»e  treated  in  f*rst  Born  approximation, 
ir verse #meari  tree  path  differential  in  momentum 
■ ' 'er,  M,  and  energy  transfer,  Ti>. , for  a particle 

• velocity  v is  qiven  by 


- j i *t  ..'esc  r i p 1 1 on  of  t 'e  nteradion  of 

e e.tTcir-  « ’h  r.atte*  over  a 'arge  ranee  of  energies  is 

• ,bec’  . * .v  ■'portance  in  a wiae  variety  of 
meme!  i a ' • d applied  areas  f rom  toe  theoretical 

tand;.  mt  »’  actions  c*  energ.  less  and  range  of 

e'ectron  rarv  different  materials  nave  formed  the 

•ai’f  * ?f  eas:  * #(  e»tef's-ve  tabulations.  ■ . Both 
t^ese  *or»  are  restricted  to  electron  energies 
keV  mj  are  v ed  n *fe  Be  the  theory  of  Stopping 
*er  'r  '.id  mg  va-  >,  "tc.‘ ' ‘ -cat  mir  arv,  corrections 

r : ge* ' t>-effe  • rre:t‘ons:  r,e  *eel  fhat  cal 

'a*iors  '"ear  ;aths,  energy  loss,  et  for  elec 

tf,’ps  of  energy  ke»,  based  or.  the  best  currently 

i va^able  theoret'm’  • oe  ‘ *• , will  provide  useful 
guides  for  interprets! ior  of  experimental  data  as  well 
as  input  for  ale ulat ions  in  applied  areas 

Edr‘.er  tabuloiicn  , • inverse  mean  free  paths, 
lopping  rowers,  ’.A  rar  :e  , and  straggling  for  Al 
and  A ;_f.  and  S'  and  c ••  C involved  tne  use  of 
.ar*ous  •iOdels  v describe  the  exte'rded  electron  states 
• a solid  valence  ban}  f-  onductioo  band',  and 
tneoreti  al  calculations  of  atone,  generalised  osc i 
'ator  ftrengths  ■ AOS ' s > to  describe  excitation  of 

• ' ec t r :,r s frorr.  the  • iner  hells  of  the  atoms  in  the 

ds  r.nf.re  deta  lea,  a.  prjo_ri  GOS  calculations  such 

in  Re1  • • iva i \ able , • tst  depend 

a . a * '■ ' based  on  classical  binary  collision  cross 

.‘•TV  r,y.ir  . *a  V<--  V .'.'!’  On  S . V approximate 

- ne  e*c i ' i * ‘ oe  the  inner -she  1 electrons.  Althouqn 
‘n.?  ’r  Is  are  not  verv  i . ,*‘tant  in  calculations 

V verse  "f^n  free  paths  f electrons  for  energies 
*ev.  the  mner-snell  electron!  do  contribute  signi 
a^t  . *•  *he  stepping  iwer  frr  incident  electron 

energies  -500  eV.,»'*'fc 

iTS’eed  of  dividing  th*  °1»c*rons  w the  solid 
r- * -1  two  dietinct  groups  t . i m ner-shell  and  valence 
band,  we  wi'l  present  here  an  alternative  approach  in 
which  the  response  of  all  of  the  electrons  in  the  solid 
are  treated  in  a ;or,prehensi ve  manner  within  the  frame 
work  of  a smqle  mode.  A model  that  has  been  devel- 
oped for  this  purpose  ?s  the  electron  qas  statistical 
model  , or  statistical  model'  for  brevity. 


•This  researv  wa  . sponsored  in  ^art  by  the  Defense 
Nuclear  Agency  under  sutifask  TA"40  and  in  part  by  the 
Energy  hpsearch  and  Development  Admtnistrat ion  under 
contrac*  with  union  Carbide  Nuclear  Corporation. 

•Postdo.  tor*  fellow,  University  of  Tennessee , Knox v f Me 
ennessee  37919 

• Computer  Science*.  Division,  Oak  Ridge  National  labora- 
tory, Oav  Ridge,  Tennessee  37830 
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Abe*e  ) is  the  exact  dielectric  response  function 

‘ the  solid  • 8 We  assume  in  this  work  that  the  solid 
isotropir  and  homogeneous 


'or  our  calculations  of  inverse  mean  free  paths 
an.:  stoppmq  powers,  it  is  sufficient  to  compute 
■nverse  mean  free  paths  differential  in  energy  transfer 
only  This  differential  inverse  mean  free  path  (DIMFP) 
for  energy  loss  by  an  electron  with  energy  F c mv2 /2 
m the  solid  is  given  bv 


were  hk.  xTr  [ *T  • *T-fL  ] and  a0  2 Tl‘  /mev  This 
•pression  assumes  that  the  energy -momentum  relation 
for  a swift  electron  in  the  solid  does  not  differ 
appreciably  from  that  of  a free  electron  In  vacuum. 

Given  c(k,u)  for  the  solid,  the  quantities  of 
interest  here  follow  directly  from  t(F,TL>).  The 
inverse  near,  free  path  of  the  electron,  u,  Is  given 
by  integrating  over  allowed  enerqy  transfers  as 


l(E)  = j d(tlw)  .(E.Tlu,).  (3) 

The  rate  of  enerqy  loss  of  the  electron,  or  the  stop- 
o »ng  power  of  the  medium,  is  given  by 


S(E)  dE/dx  = j d (Tlw ) Xu.  i (E,Xu>).  («) 

Since  exact  dielectric  response  functions  for 
solids  are  not  available,  various  models  and  approxi- 
mations must  be  employed  to  calculate  quantities  which 
may  be  compared  with  experimental  measurements.  For 
example,  the  electron  gas  model7,efor  the  response  of 
the  conduction-band  electrons  in  metals  has  been  used 
quite  widely.  The  um'prrr,  electron  density  assumed  in 
this  model  leads  to  results  in  many  cases  which  agree 
quite  well  with  experiment.  In  the  case  of  aluminum, 
for  example,  calculations  of  electron  Inelastic  mean 
free  paths  yield  excellent  aqreement  with  experiment.’ 
'he  model  which  is  used  here,  the  electron  gas  statis- 
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tlcal  model,  incorporates  the  relative  simplicity  of  the 
electron  gas  model  but  is  applicable  to  a larger  > lass 
of  meterials.  Such  a statistical  approach,  pionee  * d by 
Jens  Lindhard,  has  .met  with  considerable  success  when 
applied  to  the  calculation  of  charqed-partic le  energy 
loss  In  many  different  materials.  Basically  the 
approach  involves  the  assumption  that  the  stopping 
medium  is  characterized  by  a space- vary ing  density  of 
electrons,  n(r),  and  that  the  contribution  gf  those 
electrons  in  a small  volume  element  d‘r  at  r to  the 
stopping  process  is  the  sane  as  the  same  number  ,*  *•••♦• 
electrons  in  an  electron  gas  at  the  sane  density  '*• 
total  stopping  power, ^for  example,  is  then  compute*  by 
averaging  over  the  n(r)  distribution  appropriate  to  a 
given  medium.  The  contribution  of  the  inner -sne  ’ ele* 
trons  to  the  electron  energy  loss  is  included  in  j 
statistical  fashion  in  the  averaging  process 

The  DIMFP  for  an  electron  gas  of  umfo'1*  density  r, , 
t(E,Tiw,n),  is  calculated  from  £ } ;2'  and  used  t eval.* 

ate  u(E*n)  and  $(E,n).  These  quantities  are  • ' /ated 
for  a large ^number  of  values  of  n 'or  ^ ,e  m the  average 
over  the  n(r)  distribution.  Tne  inverse  -’ea'*  free  oath 
<-(E)>  and  stopping  power  - S( F ) * calculated  ■ "e  basis 
of  this  statistical  model  are  obta’ned  from 


-u(E) ' I d3r  v (f ,n(r) ) d . (5) 


and 

<S(E)>  djr  $(F,n(r))  , dV  (6) 

where  the  integration  is  arriel  out  over  some  appropri- 
ate volume  of  the  solid,  "he  electron  gas  model  will  be 
discussed  briefly  in  Part  B,  and  the  spherically  sym- 
metric electron  density  distributi  n(r)  calculated 
in  a Wlgner-Seitz  cell  using  a relativistic  Hartree- 
Fock  method,  which  were  used  to  evaluate  the  averages 
Indicated  in  Fqs.  (S)  and  ,6),  will  he  described  briefly 
in  Part  C. 

B Electron  Gas  of  Un  i *orm  Den  si  ly 

The  dielectric  response  function  o(  a uniform  den- 
sity electron  gas  in  the  form  derived  by  Lindhard'  fas 
been  used  extensively  The  equation  for  x,  ' an  be 
found  in  Reference  ' and  will  not  be  repeated  here 
Instead  we  present  graphical U an  example  f the  results 
of  calculations  required  for  the  averages  indicated  in 
Fqs.  (5)  and  ^6)  In  figure  1 *e  >.hOw  'he  stopp  ng 
po**r  in  .itomic  unit-,  .is  a function  ->f  int  ident  electron 
energy  for  electron  gases  ,>♦  several  uniform  densities 
Ihe  curves  are  labeled  by  the  one  electron  radius 
which  Is  relates  to  the  um for o density  n by  1/n 
4wa^rj/3  where  ao  h-  me  0.S29  lhe  stopping 

power  in  atomic  units  may  he  >. .qnverted  to  units  o'  ev  a 
by  multiplying  by  27.2  J.529  SI  4 Similar  calcula 
tlons  have  been  carried  out  for  . ’the  limits  on  the 
integrations  In  Fqs  (3)  and  (41  for  these  calculations 
are  given  by  0 <Tiw  E-Er(n),  where  Fr(n)  is  the  Fermi 
energy  for  the  uniform  electron  density  n In  terms  of 
the  'one  electron  radius  . the  Fermi  energy  is  oiveri  by 
Fr  (e2/?a0)  (9,/4)‘/*  r s • 

(.  Calculation  of  n(rj 

For  the  calculation  of  electron  tensity  as  a func- 
tion of  position  In  the  solid  we  have  used  an  exist mq 
computer  program  which  employs  relativistic  self-con- 
sistent - f lei d calculations  to  evaluate  the  spheri 

cally  Symmetric  electron  density  n(r)  in  a Wiqner 
Seitz  sphere  around  an  atomic  nucleus  As  an  example 
of  this  type  of  calculation  we  show  in  f iqure  2 the 
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Fig  I Stopping  power  of  an  electron  gas  for  several 
values  of  uniform,  electron  gas  density  as  a 
function  of  incident  electron  energy 
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' l .'  tie.  tron  densities  in  Au  and  Aq  computed  by  a 
relativistu  Hartree-Fock  procedure  usinq  thp 
wioner-Sei t:  boundary  condition. 


lensity  Anr-'n(r)  as  a function  of  r for  Ag  and  Au. 
These  densities  represent,  approximately,  the  number 
of  electrons  per  unit  radius  of  a spherical  shell  of 
radius  r centered  about  each  nucleus  in  the  respective 
solids  Atomic  units  are  used  here,  e.q.,  all  lengths 
are  expressed  in  units  of  a0  Given  these  computed 
electron  densities,  the  inverse  mean  free  path  and 
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Compar i 
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Ot he  i alculatlons 

As  jn  n i 1 1 a • test  of  the  model  described  here  we 
- pared  the  near  ‘ r»-e  piths  and  stopping  powers  calcu- 
lated * r fll  and  i «'th  >r  earlier  tabulated  re- 
sults For  tne  ean  free  paths  we  find  the  predic- 


tions of  the  statistical  model  to  be  ‘101  smaller  than 
our  earlier  more  detailed  calculations  In  the  energy 
region  100  rV  to  1 0‘*  eV  More  spec  1 f 1<  al  1 y the  elec 
r»on  mean  free  path*.  4n  St  are  11  less  than  the  values 
predicted  front  F<#»  f 4 for  f»0  eV  •{  10**  eV , while  for  Al 
we  * i r*.  1 »•)»•»  iron  mean  free  paths  Si  less  than  the 
value-,  pin!  tried  « v «wr-  U+1  1 for  100  »-V  - 1 - 10“  aV  for 

ent-r.jii-s  1 00  eV  . lar«jr»  differences  are  found  as  we 
, to  lower  energies  At  I * 10  eV , the  results  here 
ar»-  SUV  of  our  earlier  results.  1fi4  This  difference  Is 
due  principally  to  the  neolect  of  exchange  effects  be  • 
■weer  the  incident  electron  and  electrons  of  the  medium 
■>  the  statistical  model  For  the  stopping  power,  the 
statistical  model  predictions  are  ,201  larger  than  our 
earlier  results  for  100  e\  ^JO**  eV  with  larger  dif- 
ferences occur  inq  for  1-100  eV,  aqain  due  to  neglect  of 
exchange  effects.  A comparison  with  our  earlier,  pub- 
lished results  for  the  stopping  power  of  Al  not  Includ- 
ing exchange  effects1  shows  agreement  to  within  mlOX 
for  the  broad  energy  range  10  eV  to  104*  eV. 


We  show  in  Figures  3-5  the  stopping  power  of  Au , 
Aq,  Cu,  and  Ni  calculated  f-om  the  statistical  model; 
results  from  Rethe-Bloch  theory*’  are  shown  for  E>10**eV 


on**.- owe  *>  «o*4 


electron  fNC«GY  above  £f  («V) 


' ; 3 Stopping  power  of  Au  and  Aq  for  electrons  cal 

culated  using  a statistical  model. 


oRNi-owc  r*  *o<a 


Fin.  4 Stopping  power  of  Cu  for  electrons  calculated 
using  a statistical  model. 
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ig.5  Stopping  power  of  Ni  for  electrons  calculated 
using  a statistical  model. 


or  these  four  solids  the  stopping  powers  calculated  on 
he  basis  of  the  statistical  model  agree  with  the  pre- 
lctions  of  Bethe-Bloch  theory*  to  better  than  10%  at 
-10  keV;  (Ni,  9*  larger;  Cu  and  Ag,  3%  larger;  Au,  4%  , 
arger). 

for  the  energy  range  200  eV  ^t<2000  eV,  mean  free 
laths  have  been  calculated  by  Penni:  using  an  electron 
las  model  including  exchange  and  correlation  effects  in 
in  approximate  way,  and  with  approximate  corrections 
'or  inner-shell  contributions.  The  best  agreement  with 
;he  statistical  model  calculations  is  found  for  A1 
vnere  Penn's  results  are  5%  smaller  at  E<200  eV,  with 
the  difference  decreasing  to  a constant  value  of  ^3% 
for  E^IOOO  eV.  For  Si,  Penn's  mean  free  path  is  12% 
smaller  at  E=200  eV  and  '-8%  smaller  for  E>600  eV.For 
the  other  solids,  for  the  same  energy  range,  Penn's 
results  are  smaller  than  ours  by  -v20%  for  Ni,  Au,  and 
Ag,  and  ^30%  for  Cu.  We  display  in  Figures  6 and  7 
the  results  of  Our  statistical  model  calculations  for 
the  inverse  mean  free  path  of  electrons  in  Cu  and  Ni . 


OWNL-0W  7«-«» 


Fig.  6 Inverse  mean  free  path  of  electrons  in  Cu  cal- 
culated using  a statistical  model.  Experi- 
mental measurements  of  attenuation  length  were 
taken  from  the  references  shown  in  the  figure. 
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ELECTRON  ENERGY  ABOVE  £r  («V) 

Fig  7 Inverse  mean  free  path  ot  electrons  in  Ni  cal- 
culated usinq  a statistica.  model.  Experimental 
measurements  of  attenuation  length  were  taken 
from  the  references  shown  in  the  figure. 


Also  shown  are  experimental  measurements  of  attenuation 
lengths  taken  from  the  references  indicated  on  the 
figures. 


IV . Conclusions 

We  have  described  a statistical  model  which  may  be 
used  with  some  confidence  to  calculate  mean  free  paths 
and  stopping  powers  for  elemental  solids  where  more 
detailed  calculations,  as  e.g.  in  References  3 and  4, 
are  unavailable.  For  the  region  of  electron  energies 
clO"  eV,  reliable  determinations  of  Such  electron  in- 
teraction parameters  are  important  for  studies  of  low 
energy  electron  diffraction,  electron  transport,  etc. 

A more  extensive  tabulation,  including  ranges  and 
straggling,  for  this  energy  region  for  Ni,  Cu,  Ag,  and 
Au  is  available. 13 
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METRIC'.  SYSTEM 


BASE  UNITS: 

Quantity 

length 

mass 

time 

electric  c urreiit 
thermmh  nanm  temperature 
amount  of  substain  e 
luminous  intensitN 

SUPPLEMENTARY  UNITS 

plane  angle 
solid  angle 

DERIVED  UNITS: 
Acceleration 

ac  tivity  (of  a radioac  five  sourc  *•  I 

angular  a<  celeratmn 

angular  veloc  ity 

area 

density 

electric  capacitance 

electrical  c undue  tarn  *• 

electric:  field  strength 

electric  inductance 

electric:  potential  difference 

electric  resistance 

electromotive  fore  e 

energy 

entropy 

force 

frequency 

illuminanc  e 

luminance 

luminous  flux 

magnetic  field  strength 

magnetic  flux 

magnetic  flux  density 

magnetomotive  force 

power 

pressure 

quantity  ot  electric  it v 
quantity  of  heat 
radiant  intensity 
specific,  heat 
stress 

thermal  < ondur  tivity 
velocity 

viscosity,  dynamic 

viscosity.  lonematM 

voltage 

volume 

wavenumber 

work 


Unit 


metre 

kilogram 

second 

ampere 

kelvm 

mole 

i .indel.i 


radian 

steradian 


metre  per  second  squared 

disintegration  per  second 

radian  per  sec  ond  squared 

radian  per  second 

square  metre 

kilogram  per  cubic  metre 

farad 

siemens 

volt  per  metre 

henry 

volt 

ohm 

volt 

joule 

joule  per  kelvm 

newton 

hertz 

lux 

( sndela  per  square  metre 
lumen 

ampere  per  metre 

weber 

tesla 

ampere 

watt 

pascal 

c oulornb 

joule 

watt  per  steradian 
joule  per  kilogram-kelvin 
pasc  a I 

watt  per  nietre-kelvin 
metre?  per  second 
pasc  al  second 
square  metre  per  sec  ond 
volt 

c ubic,  metre 
reciprocal  metre 
joule 


SI  PREFIXES 


Multiple  ation  Factors 


j non  non  non  ooo 

HP* 

i non  ooo  ooo 

10'' 

] ooo  ooo 

io* 

1 ooo 

101 

100  - 

IIP 

10 

10’ 

0 1 

10  • 

0 01 

10  1 

0 001 

10  1 

n ooo  ooi 

10  * 

o ooo  ooo  001 

10  *' 

o non  ooo  ooo  001 

10  '* 

o ooo  ooo  ooo  ooo  001 

10  " 

non  non  ooo  ooo  ooo  nni 

10  •" 

SI  Symbol 

m 

tg 

.s 

A 

k 

mol 

cd 


rad 

sr 


\ 

S 

u 

V 

V 

I 

N 

Hz 

lx 

Im 

YVb 

T 

A 

W 

Pa 

C 

I 


Pa 


V 


) 


Formula 


m/s 

(disintegration)^ 

rad/s 

rad/s 

m 

kgm 

A-s  (X 
A/V 

V m 
V-s/A 
YV  A 

V A 
YV/A 
N-m 
I'K 

kg-m/s 

(cyc:le)/s 

lm/m 

cd/m 

cd-sr 

Aim 

Vs 

YVb/m 

)/s 

N'm 

A-s 

N-m 

W/sr 

I'kK-K 

N/m 

YV/m-k 

ms 

Pa-s 

m/s 

YV/A 

m 

( wave)  m 
N-m 


Prefix 

SI  Sym 

tera 

T 

K*Ka 

C 

M 

mega 

kilo 

k 

hecto* 

h 

deka* 

da 

dad* 

d 

c enti* 

i: 

inilli 

m 

micro 

P 

nano 

» 

pic.n 

femtn 

r 

atto 

a 

* Id  hi  avoided  where  possible 
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MISSION 

of 

Rome  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  in  command,  control,  and  coimunlcatlons 
(C3)  activities , and  in  the  C3  areas  of  information  sciences 
and  intelligence . The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability , maintainability  and 
compatibility. 
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